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MEASUREMENTS AND A MODEL FOR CONVECTIVE VELOCITIES 
IN THE TURBULENT BOUNDARY LAYER 

I. INTRODUCTION 


The extraction of an exact solution from the equations of motion for the turbulent 
boundary layer has proven to be difficult because of the lack of an adequate physical model 
for the turbulence and because of mathematical complexities. A basic step forward would 
be an understanding of the fundamental physical structure of the flow. The examination of 
organized fluid structures within the boundary layer would give valuable insight into the 
understanding of the physical flow structures. The velocities of such organized structures 
will be designated as convective velocities. 

Klebanoff [ 1 ] , Kline et ah, [2] , and Grant [3] give some insight into the turbulent 
structure. G. I. Taylor [4] with what is now called “Taylor’s Hypothesis” gave perhaps the 
first physical model for the motion of eddies embedded within a flow. Taylor stated, “If the 
velocity of the air stream which carries the eddies is very much greater than the turbulent 
velocity, one may assume that the sequence of changes in u ' at the fixed point are simply 
due to the passage of an unchanging pattern of turbulent motion over the point, i.e., one 
may assume that u'= 0(t) = 0(X/V m ) where X is measured upstream at time t = 0 from 

the fixed point where u'is measured.” Taylor requires that the turbulent structure of the 
flow is convected or transported at the same rate as the mean flow, when the above 
conditions are met. One statistical method by which to evaluate the rate at which turbulent 
structures are convected within a flow is the space-time correlation of signals originating at 
spatially separated sensors. The space-time correlation comparing the signature 
identification of turbulent structures sensed at one location with signatures sensed at a 
downstream location yields a most probable transit time for structures to be convected 
between the sensors. The rate at which turbulent structures are convected can be evaluated 
by dividing the probe separation by the most probable transit time evaluated from the 
space-time correlation. In this work the term convective velocity is defined as being the rate 
at which the turbulent structures are convected. 

The convective motion of the embedded eddies has been studied in grid-generated 
turbulence. The most recent study is that of Comte-Bellot and Corrsin [5] , 1971, where the 
grid flow was found to be consistent with the constraints of Taylor’s hypothesis. The 
experimental results were also consistent with the hypothesis. 

Fluid motions with high turbulence intensities, or those where velocity gradients 
allow the mixing of fluids of different velocities, cannot be considered to fall within the 
constraints of Taylor’s hypothesis. Thus, convective motions in shear flows, such as jet flow 
and the turbulent boundary layer, become somewhat more complex. Wills [6-9], Favre, 
Gaviglio, and Dumas [10-15] have shown that convective velocities in the turbulent 
boundary layer may be much different from the pointwise mean velocities. 

. This report presents a physical model which describes the convective velocity 
processes for the flat plate turbulent boundary layer. Analytical inferences and restrictions 
are derived from the physical model. Finally, experimental evidence is presented to justify 
the model. 
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The model predicts convective velocities greater than the mean flow near the wall of 
a turbulent boundary layer, below the formation zone. Convective velocities less than the 
mean flow are predicted for the outer portion of the layer, above the formation zone. In 
special instances the convective velocity may be greater than the mean velocity throughout 
the layer. 


II. BACKGROUND 


Osborn Reynolds [16] gives the best history of fluid mechanics before 1900. 
Reynolds notes the original accomplishments of Navierin 1822 with updating by Venant in 
1846 and Stokes in 1845, with Stokes making mention of the term “eddies.” Osborn 
Reynolds introduced what is now termed the Reynolds Number in 1883, and expanded on 
it in 1895 [ 16] . This work laid the foundation for the understanding of the general effect of 
eddies in fluid motions. In 1915, G. I. Taylor examined eddy motion in the atmosphere and 
made the first documented statement about the convective transport ability of the eddy, 
when he said, “It seems natural to suppose that eddies will transfer not only the heat and 
water vapour, but also the momentum of the layer in which they originated to the layer 
with which they mix” [ 17] . 

In 1917 Taylor performed experiments to map the shape of eddies in the 
atmosphere [ 18] . Taylor found that the eddy shape is dependent upon the elevation above 
the ground. He expanded further on the eddy in 1918 [ 19] , and in 1921, he introduced the 
idea of a random migration (walk), and with it the ability of the turbulent fluid to diffuse 
fluid properties [20] . In 1923 he introduced the concept of the eddy cascade [21 ] ; that is, 
that the large eddies lose their energy by forming smaller eddies, and at last the smallest 
eddies are dissipated by viscosity. In 1931, he examined the motion of eddies near an 
instability [22]. Taylor’s microscale and macroscale were defined in 1935 to give insight 
into eddy size as well as a statistical approach to the problem of turbulence [23-27], He 
published the relationship between correlation measurements and spectra in 1938 [4] . This 
study laid the foundation for future work using correlation measurements to describe 
particular motions of selected portions of the spectrum. It is here that Taylor gave insight 
into the motions of eddy structures as compared with mean fluid motions. Taylor states, “If 
the velocity of the air stream which carries the eddies is very much greater than the 
turbulent velocity, one may assume that the sequence of changes in u ' at the fixed point 
are simply due to the passage of an unchanging pattern of turbulent motion over the point, 
i.e., one may assume that 


u' = 0(t) = 0(X/V m ) , (7) 


where x is measured upstream at time t = 0 from the fixed point where u is measured. In 
the limit when u '/V m ->0 (7) is certainly true. Assuming that (7) is still true 

when u '/V m is small but not zero, R x is defined as 
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Ht) 0(t + X/V m ) 


This quote was later to become known as Taylor’s hypothesis. Taylor referenced work by 
Dryden [28] where correlations were measured with band filtered circuits. Dryden’s test 
was performed in grid-generated turbulence and was held within the constraints of Taylor’s 
hypothesis. Since Dryden’s experimental results which supported Taylor’s hypothesis were 
published, many others have examined convection in grid-generated turbulence [29-32] . A 
comprehensive study is reported by Comte-Bellot and Corrsin [5], and their experiments 
are shown to support Taylor’s hypothesis. 

Since it does not fall within the constraints of Taylor’s hypothesis, the turbulent 
boundary layer does not have a simple space-time relationship with which to speculate on 
the turbulent structure. Visualization and analytical techniques have been used in an 
attempt to understand the boundary layer’s physical structure. Numerous papers which 
delve into the structure of turbulence have been written. The most important studies with 
regard to the background and development of this paper are given in References 1,2,3, and 
33 through 47. 

Insight into the flow structure has been obtained by studying the structure of the 
turbulent pressure field on the floor of a boundary layer [48-52] . Corcos [48] , Bradshaw 
[49], Blake [50], Willmarth and Roos [51], and Willmarth and Wooldridge [52] note that 
the convective velocity of the pressure fluctuations is a function of wave number. Willmarth 
and Wooldridge [52] note, “The lower convection speeds are measured when the spatial 
separation of the pressure transducers is small, or when only the pressure fluctuations at 
high frequencies are correlated.” This observation also agrees with Corcos [48] . Willmarth 
also notes that the pressure fluctuations are produced by eddies convecting at the same rate, 
intimating that the convection velocities of eddies near the wall are a monotonically 
decreasing function of frequency. 

The use of correlations for resolving the structure of the turbulent boundary layer 
has not been limited to pressure measurements [10, 11]. The use of such techniques to 
establish convective properties of the flow is of prime interest for this report. Favre et al. 
[10-15] , Wills [6-9] , and others [53-58] have approached the problem of convection with 
the use of space-time correlation concepts. The constraints, limitations, and use of these 
techniques are presented in Section IV to give physical and analytical insight into a simple 
model to predict gross characteristics of the convection processes in the turbulent boundary 
layer. These concepts are also used for constructing an experimental test procedure to check 
the model’s predictions. 

Kline et al. [2] , presented extensive visual studies of the turbulent boundary layer 
near the wall, where the major portion of the turbulent production occurs. It was observed 
that well organized spatially and temporally dependent motions in the viscous sublayer 
formed low speed streaks in the region near the wall. These streaks begin a gradual “lift-up” 
process, which finally oscillates and ejects fluid into the outer portion of the boundary 
layer. Kline used a hydrogen bubble technique which showed details of the trajectory away 
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from the wall. He attributes the streak breakup to the- dominant transfer of momentum 
between the inner and outer region. These concepts appear consistent with the findings of 
Sandborn [38], who observed that the high frequency components are intermittent 
throughout the layer. . 

Laufer and Narayanan [59] gave reasons to believe that these inner streaks and their 
final breakdown may be correlated with the undulating outer edge of the flow. Kibens and 
fCovasznay [34] also noted that the streaky structure near the wall was closely associated 
with the undulating outer edge of the boundary layer. Grant [3] mentions that stress 
relieving motions very near the wall would cause an outward migration originating near the 
wall. When Kline and Lahey [54] later used short time averaging to better examine the 
migrations, the spectra of the migrations seemed to peak away from zero. One is now 
tempted to speculate about whether the migration results from fluid being driven upward 
from an instantaneous localized separation. 


III. THE PHYSICAL MODEL 


Based on the observations of Kline et al. [2], Laufer and Narayanan [59], and 
others, we may hypothesize a physical model for the turbulent boundary layer. A portion of 
the model relies on the concept that a relationship exists between the external potential 
flow and turbulent production near the wall of a turbulent boundary layer. This concept 
comes from Laufer and Narayanan [59] who showed that processes occurring near the 
viscous sublayer scaled with outer flow parameters of a boundary layer. 

The scaling may be a result of strong localized pressure gradients created in the 
valleys of the convoluting (intermittent) outer edge of the boundary layer [34, 54] . These 
strong localized pressure gradients in the outer portion may create either small jets or 
pressure pulses which travel through the boundary layer and result in an instability near the 
surface. Directly above the viscous subregion near the wall is the layer where the effect is 
observed. The subsequent instability becomes the streaks noticed by Kline which form 
packets of vorticity. The packet of vorticity, with its circulation, is then acted upon by the 
surrounding velocity field. The direction of circulation coupled with the surrounding fluid 
velocity then creates a pressure gradient across the packet, causing the packet to be driven 
either upward or downward toward the wall. If the surrounding fluid’s velocity coupled 
with the circulation velocity on the top is greater than it is on the bottom, the packet will 
be forced upward. Conversely, if the surrounding fluid’s velocity coupled with the 
circulation velocity on the bottom of the packet is greater than the top velocity, the packet 
will move downward toward the wall. 

A model allowing turbulence to enter the viscous subregion tends to be in agreement 
with Klebanoff [ 1 ] , who notes that the dissipation is greater than the production in a very 
thin region next to the wall. (The above model allows the packets of vorticity to enter the 
viscous sublayer where they may be quickly dissipated.) After the packets of vorticity have 
been formed, they may be identified by the structure of their oscillations. The velocity of 
the packets is termed “convective velocity” because the mean flow becomes the forcing or 
driving function on the packets and convects the packets in a downstream direction. 
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Corcos [48] , Bradshaw [49] , Blake [50] , Willmarth and Roos [51], Willmarth and 
Wooldridge [52] and Wills [6-9] all observe. that the convection processes near the wall are 
faster than the mean motion. Experimentation by Kline et al. [2, 33, 54] shows mainly the 
motion away from the sublayer. Favre et al. [10-15] found that the convective motion in 
the outer portions of the layer is slower than the mean motion. 

The input of energy from the highly intermittent outer edge of the boundary layer 
will be a rapid process. The packet of vorticity is formed and moves approximately with the 
mean velocity of the layer in which it is formed. The packet of vorticity moves both upward 
into the main body of the boundary layer and downward into the viscous subregion. The 
region directly above the viscous subregion is where Kline [2] and Kovasznay [42] both 
observe the apparent formation of the vorticity. 

The present experimental study will show that the convective velocity is equal to the 
mean (ensemble average) velocity in the region where the packets are formed. Below the 
formation region (viscous sublayer), the convective velocity is greater than the mean 
velocity; and above the formation region, the convective velocity is less than or equal to the 
mean velocity. (The relative positions of the various regions are shown in Figure 1 .) 

1 FREE STREAM 



The packet of vorticity, now an integral part of the boundary layer, is traceable 
because of the packet’s organized wave pattern. The high degree of orderliness of the 
structure makes it identifiable for distances of two orders of magnitude greater than its 
length scales [60] . This characteristic suggests that once the packets are formed they tend 
to be entities within themselves. The present mathematical analysis of the proposed model 
will show that the intermittent outer edge of the boundary layer is controlled by the 
packets, which are shown to control the asymptotic growth of the boundary layer. Insight 
into the question of when a flow is fully established is also gained through the present 
mathematical study. 
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The derivations for the equations of motion of the packets of vorticity will be 
presented later, and will be inspected to obtain inferences regarding the motion of the 
packet. 


A. Vertical Direction of Motion 

The following assumptions are made so that an analytical derivation can be 
performed: 

1. Once the packet leaves the layer in which it was formed only its translational 
motion heed be considered. 

2. Since the analysis starts only after migration from the production zone,' it may 
be assumed that the packet has an instantaneous velocity at time t = 0 . 

3. Once the packet has moved away from the production zone the only 
accelerational influence on the packet in the vertical direction will be opposite to its 
velocity, such as a pressure drag or viscous drag. 

Consider a packet of vorticity that is three dimensional. Three characteristic lengths 
or scales, Xj, X 2 and X 3 , are used to describe the packet. The mass of the packet is 
represented by C 1 X 1 X 2 X 3 p. (C, = tt/ 4, if the packet were considered a cylinder.) 
Subscripts 1, 2, and 3 denote directions parallel to the mean flow, vertical to the mean flow 
away from the wall, and parallel to the wall and perpendicular to the mean flow, 
respectively. 

When the previous assumptions for the physics of the packet’s vertical motion are 
used, the only force acting on the packet will be in a direction opposite its velocity. The 
effect of drag is introduced by using a viscous damping term, p , with the dimensions 
FT/L 2 . If the flow were completely viscous, p would be some constant times the dynamic 
viscosity. In dimensionally forming Newton’s second law, a characteristic length is 
necessary. The characteristic length perpendicular to the direction of motion was selected 
for the following analysis. However, the conclusions derived from the following analysis do 
not depend upon which characteristic length is selected. This is because the analysis 
considers all packets to have similar geometrical configurations; thus each characteristic 
length is directly proportional to the other characteristic lengths. The reason we use Xj for 
the characteristic length is because Cliff [60] used this characteristic length to formulate an 
empirical formula describing radial motion of turbulent packets in an axisymmetrical jet. 
The motion was successfully used to predict the intermittency factor as a function of radial 
direction in the axisymmetrical jet. Noting that the acceleration is in a direction opposite to 
that of the velocity, we see that Newton’s second law becomes 


Ci pX i X 2 X 3 


d 2 y(t) 
dt 2 


+ pX, 


dy(t) 

dt 


= 0 


( 1 ) 


6 



or 


d* y(t) + _ P dy(t) _ n 

dt 2 C,pX 2 X 3 dt u • (2) 


The physical restrictions placed on the model require that at time, t , equal to zero, 
an initial velocity exists, Vj, and the original location may be some vertical distance above 

the flat plate; i.e., y(0) = Yj. A reasonable assumption for the region of formation might 
be y* = 50 [61], where y*=yu 7 ./i?. Here, u T = shear velocity = V r wa j j/p and 

t \ = kinematic viscosity. The initial conditions may be expressed as follows: 


dy(t)/dt = V y j 


and 


(3) 


y(0) = Yj , 


at t = 0. 

Integrating equation (1) and using y(0) = Yj yields 


dy(t) 

dt 


P 

C,pX 2 X 3 


y(t) + Vyj 


+ 


5 Op 2 
Ci x 2 x 3 


^w P 


(4) 


Using the translated coordinate system with y = 0 at y* = 50, equation (4) becomes 


dy(t) 

dt 


P 

Ci pX 2 x 3 


y(t) + Vyj 


(5) 


Equation (5) is now the equation for the velocity of the packet of vorticity in the vertical 
direction, the left-hand side of the equation being the vertical velocity of the packet. The 
first term on the right is the dissipative term which retards the motion of the packet. 
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A direct solution of equation (2) with y(0) = 0, dy(0)/dt = Vyj , at t = 0 gives 


y(t) 


y Ci pX 2 X 3 

yi p 


e 


-f—U— 

j C,pX 2 


x 3 


dt 


+ 


CipX 2 X 3 

yi p 


( 6 ) 


or 



y(t) 


y Ci pX 2 X 3 I . 

yi p ' 



p 

Ci pX 2 x 3 



(7) 


Because 


yi 


. f — 

e CipX 2 X 3 


dt 


j*y(t) 
Ci pX 2 x 3 


+ v yi 


( 8 ) 


differentiation of equation (6) produces equation (5). 


B. X-Direction of Motion 

Once a packet of vorticity has been formed and is injected into the surrounding 
boundary layer, the vertical position may be described by equation (7). 

The horizontal motion of the packet may be examined by assuming that the only 
external force in the X-direction is caused by the mean motion' of the surrounding fluid. 
That is, whenever the packet velocity is different from the mean velocity, the mean motion 
becomes the driving force to either accelerate or decelerate the packet in question. 

A first approximation might be that the force on the packet would be equal to the 
rate of change of momentum incurred by the mean fluid because of the presence of the 
packet of vorticity. The mass of fluid that must be deflected in some manner as it passes 
oyer the packet is Gj X 2 X 3 p(V m - V x )dt. A moving reference frame is used for convenience 


The force on the packet must be equal to the mass times the acceleration of the 
surrounding fluid as a result of the packet. This force must in turn accelerate the packet of 
vorticity, the total force on the fluid being • . 
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F x = C^AspCV^ - V x )dt-^3-- , 


where V x = convective velocity of packet and V m = velocity of mean motion. Since the 

actual acceleration may not be known because of the uncertainty of the shape and other 
related processes (i.e., pressure drag, etc.), a coefficient of drag is introduced such that 


C d< V m- V x> = dV m 


Thus, setting F x equal to mass times acceleration yields 


CdC.AjXaPfVjn - V x ) 2 = 


dV Y 

C1PX1X2X3 “j^— 


or 


£Xx _ C d< V m - V x ) 2 
dt X i 



2V V 
z v m v x 


V 2 1 
m J 


(9) 


( 10 ) 


In terms of displacement in the X-direction: 


d 2 x _ £d / dxV 9V £d dx 

dt 2 " X, \ dt / m x, dt 


CjV 2 
^d v m 

X, 


0 


(ID 


Solutions of these differential equations are given in Appendix A (63, 64], when the 
functional forms for C^, Xj and V m are known. When C^, X 1 , and V m are essentially 

constant, the exact solution to equations (10) and (11) are given (Appendix A). The 
solutions for constant C^, X x , V m are 


v x (t) 


*i(V m -Vi> 

m ~ tC d( V m- V i) + 


( 12 ) 
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where Vj is the initial value of V x at time t = 0 and 


x(t) = V m t - ^ Sn[Xj + tC d (V m - V^] + (13) 


The assumption that the rate of change of may be considered small can be 
justified by Kline’s [2, Fig. 15] picture of the lift-off phenomena neax the wall. Cliff [60] 
presents pictures showing shape changes in the outer portion of a free shear flow. Both of 
the above references indicate that the length scale in the X-direction may be considered 
constant for short periods of time. The rate of change of C d for cylinders and spheres is 

not greatly affected by small changes in Reynolds number, as shown by Schlichting [65] 
and Sabersky and Acosta [62]. It is shown that C d may vary approximately 30 percent 

over two orders of Reynolds number. The rate of change of V m may be considered small 

in the outer portions of the boundary layer, or where there is little shear. One could also 
develop the equations and use small increments of time to allow the equations to be used in 
high shear regions; that is, assume that the velocity takes step increases, and use the 
equations for each step. 

Once the packet of vorticity is removed from the region of high shear, the above 
simplifying assumptions are reasonable. 

The equation for the Y-direction of motion of the packet indicates that the packet 
moves most quickly away from its origin in the high shear region, and will spend most of its 
time in the outer portions of the boundary layer. Thus the above solutions should give good 
first approximations, and are exact solutions for the case where the simplifying assumptions 
are met. The limiting case where the packet is in the extreme outer portion of the boundary 
layer and is not subject to great shears should give the best results. 


C. Evaluation of the Motion of the Packet of Vorticity 

The motion of the packet of vorticity gives direct insight into the structure of the 
turbulence and the turbulent bounday layer. First, some of the packets will eventually 
become the turbulence at the intermittent outer edge of the boundary layer. Second, the 
boundary layer growth is eventually controlled by the motion of the packets, the limiting 
case occurring when the packets force the boundary layer to grow as a linear function of the 
X-direction. Third, the turbulent flow may not be considered in equilibrium as long as there 
are accelerations of the packets within the flow; however, the following analysis shows that 
• within a finite distance the packets become very close to equilibrium with the surrounding 
flow. 


The equations to be examined are as follows. 
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The acceleration of the packet in the Y-direction; 


Ay(t) 


11 V v (t) = - * . V,,; e 


c,px 2 x 3 y w c,px 2 x 3 yi 


' f C,pX 2 X 3 dt 


( 14 ) 


MV v 


(Ci 


a: _ vft A _ Jim 

ipX 2 X 3 ) 2 CipX 2 X 


The Y-direction of velocity; 


V y (t) = V yi e 


-/ 


Ci px 2 x 3 


dt 


= V. • - 


My(t) 


y 1 . Ci px 2 x 3 


(15) 


The Y-direction of displacement; 


y(t) 


Ci pX 2 X 3 

yi p 


e 



C i p X 2 X 


- dt 

3 


) 


(16) 


The X-direction of acceleration; 


A m C dt V m - V x(t )] 2 ^(V^V ,) 2 

A * w X, [tCd^-vp + X,] : 


C d^( V m- V i > 2 




| x(t) + 8n[X , + tC d (V m - Vj)] - y^(V m -V,) 


+ Xi 


(17) 


The X-direction of velocity; 
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v x (t) - V m 


X 1 (V m -V i ) 
tCd^-Vj) + X. 


= V 


m 




C d< V m - V i> 



®n[X, + tC d (V m - Vj) ] 



(18) 


The X-direction of displacement; 


x(t) - V m t -^finfX, + tC d (V m -Vj)] + - - ^ nXl 


(19) 


i 

The equations of motion for the packet of vorticity in the Y-direction, equations 
(14), (1 5), and (16), will be discussed first. 

Equation (14) requires that the deceleration of the packet be inversely proportional 
to X 2 X 3 . This inverse ratio means that larger packets are decelerated at a slower rate than 
the smaller packets. Consider the limiting case: 


limit^oo A y (t) = 0 


The acceleration becomes equal to zero when y(t) = Vyj . - 1 — — . A longer distance is 
required for the larger packets to reach a constant velocity. 

Equation (15) shows that the velocity of the packet decreases exponentially with 
time. It also requires that the vertical velocity of the packet become equal to zero at a 
C X X 

distance y = Vyj — - — J - . Thus, the penetration of the packet into the flow is 

proportional to X 2 X 3 . For the same initial vertical velocity, a larger packet will move 
farther into the flow than a small packet. Similar results were obtained by Cliff [60] for jet 
flow. One may infer from this that the spectral content of the flow should change such that 
there must be more high wave number structures near the wall where the packet forms than 
farther away. That is, the spectral roll-off must be somewhat less close to the formation 
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zone, which Tieleman [66] has shown using the CSU-Army wind tunnel. Since the 
equations are not limited to the positive y axis,- the same effect should occur below the 
formation zone as well. Thus, as one moves away from the production zone, the high 
frequency must disappear. 

Equation (16) is the general relation for the vertical displacement of the packet. At 
time t = 0, the packet is at the axis of the coordinate system; i.e., t = 0 and y = 0 where 
y = 0 is the formation zone, y* = 50 . The most interesting case is the asymptotic values for 
t very large, 


y(~) - V yi £i£hh . . (20) 


This equation implies that, if the outer edge of the boundary layer is formed by 
packets generated near the wall, only the largest structure would have propagated to this 
area. If Vyj' p , , p are . the same for two different packets, the final distance in the 

Y-direction to which they may propagate is proportional to X 2 X 3 . If X 2 and X 3 are each 
proportional to some basic characteristic length, X, the migration distance the packet will 
attain is proportional to X 2 . The larger packets are decelerated at a lower rate than the 
smaller packets. Thus, the outer edge of the boundary layer will be composed primarily of 
the larger structures. Photographs, as well as experimental results; indicate that this is true 
[33, 34, 37, 60] . The height of the boundary layer, 6 . becomes the value of y(°°) or 
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Ci pX 2 x 3 
u 


f * . ' ■ ■ . ! 

■■■ ( 21 ) 


The clumping of smaller packets could increase their 'basic size to allow the packet 
to penetrate the flow to a greater distance. Fiedler and Head [37] suggest the possibility 
that the transition zone is a breakdown of small randomly distributed patches of turbulence 
which finally coalesce. This is found to be consistent with the proposed model. For 
boundary layers, an increase in the energy of' the power spectrum at lower wave numbers 
with increasing distance from the wall could indicate the occurrence of clumping, However, 
other phenomena such as dissipation and slower migration rates of smaller packets will also 
cause a decrease in the power spectrum at large wave numbers. Tieleman [66] ; indicates a 
gain’in power in the low frequencies with increasing distance from the wall. 

Equation (15) should also be valid in the transition region. This suggests that the 
transition region must develop much more rapidly than the turbulent boundary layer. This 
may be deduced from the fact that the packets move out more quickly from the region of 
their point of formation. , . 
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Cliff [60] noted for jet flow that, since the equation for the vertical motion is 
exponential in nature, the intermittency of the outer edge should also have an exponential 
form. The present analysis for the boundary layer will lead to the exponential form for the 
variation of intermittency. Fiedler and Head [37] obtained an exponential variation of the 
outer edge intermittence for a number of turbulent boundary layers. 

The limiting case for large time gives 


limit'j'_ > oo y(t) = 5 « X 2 X 3 


( 22 ) 


(Note: limit j^. QO y(t) + distance from the wall to the formation zone = 5 . However 
limity^. 00 y(t) > distance from the wall to the formation zone.) The boundary layer 

thickness will be proportional to X 2 X 3 . Taylor [20] found that the growth of a length 
scale will be proportional to the square root of time. Thus, 


limitj^^ y(t) = 5 t^ 2 1^ 2 <* t 


Therefore, for a constant V,*, 




6 « x 


(23) 


Thus, for the limiting case of long times, the boundary layer growth becomes a linear 
function of distance. Zoric [67] points out that this is a requirement for similarity solutions 
to exist in the turbulent boundary layer. That is, one of the requirements for the boundary 
layer to have a nondimensional velocity profile that is invariant with respect to downstream 
distance is that the boundary layer have a linear growth in the streamwise direction. 

Equation (17) gives the acceleration of the packet in the X-direction. At time t = 0 , 
the acceleration is inversely proportional to its characteristic wavelength in the mean flow 
direction. As time becomes larger, such that tC^(V m -Vp >X, , the acceleration of the 

packet becomes proportional to X t and inversely proportional to time and . Thus' the 

smaller the packets are, the more quickly they are initially accelerated; that is, the smaller 
packets are the most affected by the initial conditions. Since the acceleration is also a 
function of the square of the difference in velocity from that of the mean flow, there comes 
a point where the smaller packets are moving faster than the larger ones. Thus at some time 
the larger packets (moving more slowly) are accelerated at a greater rate than the smaller 
ones. This point occurs when the finer structure has approached the mean velocity and the 
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velocity of the larger structure is still significantly different from the mean. The velocity 
equation for the X-direction is equation (18), which shows that, at time t = 0, the 
convective velocity is not a function of wavelength or characteristic length. The zero time 
limit corresponds to the formation zone. As time becomes a factor, the velocity of the 
packet of vorticity in the X-direction is negatively proportional to Xi , the characteristic 
wavelength in the mean flow direction; i.e., 


tC d( V m - V i> ► 


i 




(24) 


Thus, the smaller the packet is, the closer its velocity is to the mean velocity. Likewise, the 
larger the packet is, the greater the difference its velocity is to the mean velocity. 

Equation (24) also indicates that for t very large, the convective velocity of the 
packet once again becomes independent of wave number. At some finite time, velocities of 
all packets, which had been generated for that length of time, should travel at 
approximately the same rate. 

A simple calculation can be made for the time needed for the convective velocities 
to approach their asymptotic value. Consider that the fluid impinging on the packet is 
deflected entirely in the Y-direction. Here, then, C d = 1 . Also, when the convective 

V 

velocity is approximately equal to the mean velocity, X, would be equal to or the 

inverse of the wave number. The flow will be considered stable (i.e., a similarity solution for 
the velocity profile will exist) when the convective velocity is equal to 0.99 V m , 
where V m is the mean velocity of the outer flow. 

Equation (24) may be written as 


1 - 


v x _ X, 
V m " tV m C d 


For the outer regions, 


X. = 


Zhl 

2iri 


(25) 
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thus, 


1 



1 

27rftC d 


Applying the above requirement on the convective velocity (0.99 V m ), 


0.01 = 


2wft 


Thus the time for flow stability becomes 


t = 


1 

0.063f 


for 


V m- V x 

V m 


= 0.01 


(26) 


the general relation being 

V /(v - V ) 

t = m M m w x J 

2rrf 


(27) 


Example 1 : Consider that the freestream velocity, is = 20 ft/sec 
(6.1 m/sec) and that the initial velocity of the packet is Vj = 5 ft/sec (1 .525 m/sec). 
Consider also that the minimum frequency content of concern is 10 Hz. Therefore, 




_ 20 
"63 


0.318 ft (0.097 m) 


and 


t 


sec 


1.59 sec. 
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Equation (19) gives the distance the packet must travel before a similarity solution 
could be expected : 

x(t) = V m t - ^-entX, + CjCVhj- Vj)] + Xi 

x = (20)(1.59) - 0.318 Qn[0.318 + 1.59(15)] + 0.318(-1.15) 

= 31.8-3.2 -0.4 = 28.2 ft (8.6 m) 


For the same initial conditions, for a frequency content of 20 Hz, we see that only 
approximately one-half of this distance would be required. Equation (27) shows that the 
time required for the convective velocity to approach the mean velocity is a linear function 
of the percent difference that may be allowed. 

Example 2: 

V^, = 80 ft/sec (24.4 m/sec) , 

Vj = 1 /4 = 20 ft/sec (6.1 m/sec), 

Minimum frequency of concern = 1 0 Hz, 

Experimental error ~ 4 percent (agreement between mean and convective 
velocity), > . 

= 1.273 ft (0.338 m). ' 


, . _ 80 
and \ 2rrf 20n 


Thus, equation (25) becomes 


1 - 


V 


= 0.04 = 


m 


1 

2nftC r 


0.1 

27ft 


t = 


0.1 

(0.04)277 


= 0.398 sec 
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and 


x(t) = v m (t) + tc d (V m -v i )] 

= (80)(0.398) - 1.273 «n[ 1.273 + (0.398)(60)] + 1.273 fin 1.273 
= 27.4 ft (8.35 m) . 


This result agrees well with Zoric’s work taken in the CSU-Army wind tunnel [67]. Zoric 
finds that for V M = 80 ft/sec (24.4 m/sec), the boundary layer thickness, 6 , displacement 

thickness, S * , and momentum thickness, 6 , all begin to grow linearly at x = 30 ft 
(9.15 m). Our model is presently concerned with the linear growth of the boundary layer 
thickness. That is, when the internal flow structure comes into equilibrium with the mean 
flow, the growth of the boundary layer is predicted to be linear. 

Equation (19) gives the X-direction of motion for the packet. The larger particles are 
found to travel the least distance because their acceleration is originally somewhat slower 
than that of a smaller particle. Examples 1 and 2 show the importance of this equation. 
With the knowledge of the difference between the convective velocity and mean velocity, 
one can estimate the length required for similarity solutions to exist. The mean velocity 
profile used in the derivations may be slightly modified; however, since the nature of the 
vertical direction of motion demands that the packet will expend most of its life in the 
outer portion of the flow, equation (19) should be useable with little or no modification. 

A more detailed approach to evaluating equation (10) would be to assume a specific 
velocity distribution through the boundary layer rather than a characteristic velocity, V m . 

As an example, consider V m = [C y y + Vj] , where C y is a constant and Vj is the velocity 
in the X-direction at y = 0. 

Using a first order approximation for y(t) (i.e., V m = C y Vyjt + Vj), equation 
( 1 0) becomes 


dt 



- 2 Tf [C vV + v iJ v x + T7 Ic v v yi t + v i^ 


yi 


Xi 


yi 


(28) 


Equation (28) is of the same form as that considered in Appendix A. The method of 
approach given for the general case in Appendix A is, however, not well suited for this 
particular differential equation. Thus, a slightly different transformation is used, as outlined 
in Appendix B, to solve equation (28). 
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The solution of equation (28) in Appendix B, gives the following X-direction 
relation for the motion of the vorticity packet. A flow of the form V m = C y y , 

where C v is some constant, is assumed. The velocity in the X-direction is given by 


V x = 


-C C • 
Syi 


(Cd C v V yiAl ),/j 


tanh 


C d C v V vi 


, V2 


t + 


C V t + V- 
S yr v i 



(29) 


while the motion in the X-direction is given by 


x(t) = C y V yi -j + Vjt - ■ £n cosh 


^dCyVyiV 72 
V X > / 


V m + V i , X, f , 

~ t - ■— fin cosh 

1 L d 


*7< V m - Vf)t 


(30) 


Equations (29) and (30) are approximations to the governing equations for the 
conditions when the packet is near the floor. The velocity profile of the turbulent boundary 
layer might be examined in two parts, A and B. A is the linear fit to the portion near the 
formation zone, and B is the linear fit to the outer region. Equations (17), (18), (19), (28), 
(29), and (30) could be used to describe the motion of the packet throughout the layer. 
The asymptotic solutions would still be governed by equations (17), (18), and (19) which 
have previously been discussed. 
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A general form V m = (C y Y) ^ [62, page 229], where C v is a dimensional 

constant, may also be used for the velocity profile. Using again a first order approximation 
for the Y-direction of motion; i.e., y(t) = V y jt , gives V m = (C y V yi t) 1//7 . 

The differential equation governing the motion is still equation (10). The solution 
given in Appendix C is 


V x = 


m 


if (V m- V i } 


V m( V m - Vj) 

7 8/ 

~ C V t 7 + V t 
§'-'v v yi l v r 


(31) 


+ V 


m 


where 


v m . = (C y V yi t) ,/7 + ^ = (C v y) 1/7 + Vj 

and 


a = 2 


Equation (31) would be difficult to work with, since the velocity profile increases 
to infinity as time goes to infinity. However, the expansion of the velocity profile, 

V = Cy , in a Taylor series yields 


V(y) = V(y c ) + (y -y c )yy 6/7 - (y -y c ) 2 Cy * 3/ i 7 + . . . 


Thus in regions where y is large, the velocity profile may. be considered constant for large 
increments of elevation, implying that the equations for convective velocity in a flow with 
no velocity gradient could be used. The simple case of allowing the mean velocity to 
approach some asymptotic value and using the equations for a constant mean motion would 
be realistic. 
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* The velocity of the packets in the formation region should be independent of wave 
number, because at this location all packets move with the local mean velocity of the fluid. 
Thus, a direct transformation .between wave number and size may be determined for this 
region. When the velocity of the packets approaches the mean velocity, once again a unique 
relationship will exist between wave number and size. 

From the previous analysis one may infer that the boundary layer must grow rapidly 
through the laminar to the turbulent transition region. The boundary layer growth rate must 
then decrease with the limiting case being when the growth rate has decreased to a linear 
function of X. The shape of the undulating outer edge may be examined in the light of the 
model. The outer edge will contain, primarily, large scale structures. If the boundary layer 
goes to separation, the formation zone would be raised suddenly, and the outer edge would 
contain a greater number of smaller scale structures. One might expect that the outer edge 
would, thus, be somewhat less intermittent than upstream. (This was observed by Fiedler 
and Head [37].) 

At the point in time where the convective processes travel at nearly the mean flow 
rate, correlation measurements should more nearly indicate a near frozen flow. The width of 
the correlation curve should not change greatly with separation distance, because the 
spreading of the covariance with distance is partially attributed to the» fact that various 
frequencies traveling at different rates will not allow all frequencies to be in phase at 
separations other than zero. It may be shown that spatial weighting of individual harmonies 
will not decrease the magnitude of the normalized cross covariance (Appendix D, paragraph 
E). This leads to an interesting speculation that the majority of spreading or diffusion 
processes in a boundary layer may be caused by varying convective velocities of structures 
or frequencies of various sizes. 


IV. EVALUATION OF CONVECTIVE VELOCITIES 


This section concerns the use of the statistical ' concept of space-time covariance 
computation to obtain the magnitude of the convective velocity. The convective velocity is 
most easily obtained by either holding time constant and varying the separation distance 
until the cross-covariance becomes a maximum, or holding the separation distance constant 
and varying the time delay until the cross-covariance becomes a maximum. In both cases, 
the convective velocity is obtained by dividing the appropriate separation distance by the 
appropriate time delay. The special case of' Taylor’s frozen flow theory is covered, along 
with the case in which the flow may not be considered frozen. 


A. Convective Velocity for Time Invariant Covariance 

The relationship between time and space in shear flows has, for the most part, been 
avoided in the past by assuming that Taylor’s hypothesis [4] was valid. That is, for flows in 
which the turbulent intensity is low, the flow pattern is assumed frozen. Thus, an observer 
moving with the mean motion of the flow would detect an unchanging spatial pattern in the 
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flow. A shift in space could be represented solely by an appropriate time displacement. A 
function of space and time, f(x,t) , would be equally represented by f(x + Ur , t + r) . 
For this special case, the correlation or covariance function would be invariant in time for a 
statistically stationary set of data: that is, 


-|r f(x,t) f(x + Ur, t + r) = 0 . (32) 

01 

This result follows the analysis given in Appendix D, which requires that 


f(x,t)f(x + Ur,t + r) = f(x,t) f(x,t) = 0 . (33) 


The covariance may also be expressed as 


f(x,t) f(x + ?, t + r) = f(x,t) f(x + £ + Ur, , t +r + r, ) 


(34) 


Wills [49] noted that temporal invariance of the covariance may be obtained if it is 
assumed that the change in the turbulence pattern between the two measuring points is 
statistically independent of the instantaneous velocity fluctuation at the initial point. Wills 
obtained the relation 


d n f(x,t) 3 m [f(x,t) - f(x + |, t + r)] _ 0 

3t n 3t m 


(35) 


for some £/r and for all N and M . For n = m = 0 , 


f(x,t) 2 = f(x,t)f(x + £, t + r) 


(36) 


Frozen flow requires that f(x + |, t + r) reduce to f(x,t) when considered at the point x + % 
and at time t + r , where £ is in the direction of the flow. Thus, for frozen flow, the 
covariance is invariant. For the frozen flow, the convective velocity may be uniquely 
defined by 
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Convective Velocity. 


(37) 


^ 7 / C xy (X,t) = constant 

This relation is shown schematically in 
velocity identical to the mean flow rate. 


C v = 

Figure 2. Here the spatial pattern is convecting at a 


B. Convective Velocity for Noninvariant Covariance 

The real world presents a somewhat more complicated picture in that the flow is not 
frozen, and the covariance is not invariant. Turbulent decay, for example, is a condition that 
will cause the turbulent velocity to change its nature and the covariance to decrease in the 
streamwise direction. However, turbulent decay may or may not change the value of the 
normalized covariance (Appendix D, paragraph F). 

Figure 3 shows how the isocovariance curves may appear in nature. Dumas, Favre, 
and Gaviglio [14] present similar curves for velocity fluctuations in a flat-plate turbulent 
boundary layer. 

Two separate convective velocities may now be obtained. The first is defined by a 
two-point measurement separated by a given distance, £, where the time delay between 
sensors (located at each point) is allowed to vary. The convective velocity is equal to the 
separation distance divided by the time delay which gives the largest covariance amplitude. 

This convective velocity is represented in Figure 3 by holding £ constant at £, , and 
varying t until C v „(£,r) becomes a maximum; that is at the point 

Ay 


TT c xy«-'> = 0 


which is the case r 2 in Figure 3. Thus, the convective velocity becomes C VT = £i/t 2 , 
where C yT is the convective velocity obtained by holding % constant and varying r. 

A second convective velocity may be obtained in a similar manner by holding 
t constant and varying % . This is represented in Figure 3 by holding t 2 constant and 
moving in the % direction until C X y(?,r) becomes a maximum; that is, 


dC xy ($,r) 

9 | 


0 
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Figure 2. Iso-correlation lines for frozen flow. 


(SPACE ) 
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Figure 3. Iso-correlation lines for nonfrozen flow. 


For this case, Figure 3 gives C y £ = £ 2 /r 2 , where C y | is the convective velocity defined by 
holding t constant and varying £. 


Line A. (Fig. 3) passes through all points of 


3C xy a,r) 


3r 


= 0 . 


Therefore, the slope of line A is C VT . Line B passes through all points where * 


3C xy (| ,r) 



Therefore, the slope of Line B is C y £ . Figure 3 shows that the two convective velocities are 
different for nonfrozen flows. 

To demonstrate that the convective velocity measured using a constant spatial 
separation and varying the time delay is the more useful for the present analysis, a signal 
representative of a particular packet of fluid at several locations was considered (Fig. 4). 
Another aid is the examination of the interactions of the two separate correlations. Assume 
that the time-varying covariance between points 1 and 3 of Figure 4 is that of Figure 5. 
Let % = |i, of Figure 2, where Figure 2 gives the maximum of the curve occurring at time 
delay, r 2 . For this case, C xy (|,r) = 0.7 . From Figure 2 we see that, for the same time 

delay, we are able to get a higher correlation, approximately 0.74, if the separation distance 
was only £ 2 , (? 2 <?i ). This point is represented on Figure 5 by point A. Examining Figure 
2 again, it is seen that for £ 2 , an even larger correlation occurs at r 3 ,(t 3 <t 2 ) . This point is 
plotted as point B on Figure 5. 

Applying the same procedure we see that for r 3 there occurs a greater correlation 
at an even shorter separation, £ 3 . This correlation is point C on Figure 5. For | 3 there is a 
greater correlation at t 4 , point D of Figure 5. 

This process may be repeated over and over again until the origin is reached. Thus, 
the peak of a time covariance must lie within the time covariance curve for a slightly smaller 
separation distance. The maximum time covariance value for a given time delay (note that 
this is not the peak) is the true maximum of the space covariance for that time delay and 
corresponds to the separation distance for the curve on which it lies. For Figure 4, the time 
covariance of points 1 and 3 would give a maximum at some delay; e.g., tj . At this same 
time delay, the covariance of points 1 and 2 could have a larger value. 
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Figure 4. Schematic for correlation analysis. 


Figure 6 may be examined in the same manner as Figure 5. For a given time 
delay, r 2 , a peak will occur at £ 2 (point A, Fig. 6). However, at £ 2 a higher covariance is 
obtained at t 3 , r 3 <r 2 (point B). At t 3 a larger covariance is found at £ 3 (point C). 
At £ 3 the maximum covariance occurs at r 4 , r 4 <r 3 (point D) and so on as in the above 
case. 


Statistically, the peak of the space-time covariance curve with varying time delay 
occurs when the signals are most closely matched; that is, the optimum time delay is the 
time delay where the wave shapes are the most alike. The space-varying, space-time 
covariance curve indicates that a slight decrease in separation will give a higher covariance 
amplitude for the same time delay. The difference between the two convections gives a 
measure of the degree to which the flow is not frozen. 

The present analysis deals primarily with the passage of fluid packets and the signal, 
produced as a result of them. A transit time needed for the packet to best reproduce itself is 
desired. Thus, the time-varying scale was selected. That is, the convective velocity will be 
obtained as the given separation distance divided by the time delay at which the covariance 
becomes a maximum. (This is consistent with the work of Favre, Gaviglio, and Dumas.) 
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A typical covariance curve is shown in Figure 7. Unless otherwise stated, the covariance will 
henceforth be written as C xy (r) . The above chosen convective velocity, C VT , will be 

written V c for the overall convective velocity. To this point, nothing has been said about 
whether the convective velocity is a function of wave number or characteristic lengths. 


27 




28 


Figure 5. Correlation in time domain as function of space. 


FOR EACH CORRELATION TAKEN AT A GIVEN TIME 
DELAY, Tj, THERE IS A MAXIMUM AT SOME 


SEPARATION, ^ . AT THE SAME SEPARATION, THE 
CORRELATION HAS A GREATER VALUE AT A TIME 



Figure 6. Correlation in space domain as function of time. 


The convective velocity of specific frequencies will be denoted by Vf , where the 
subscript f will be replaced by the frequency in question. 

C. Evaluation of Convective Velocities 

The following paragraphs on concepts and experimental procedures deal primarily 
with correlation, spectra, and statistical signal analyses. A short outline on signal analysis 
and correlation concepts is given in Appendix D, which includes the spectral relationships 
for the correlations discussed in the preceding section, using both stationary and 
nonstationary series. Also included are restrictions imposed on the correlation and spectra 
by physical limitations, such as records of finite time and stationarity of the original signals 
through the nth and mth moments, respectively. The effect of correlating derivatives, a 
physical meaning for the cross phase spectrum, effects of space and frequency weighting 
functions on covariance computations, as well as the interrelationship of correlation and 
probability density, are examined. Although the fluid model presented in this report does 
not rely on these statistical tools, advances have been made which give more insight into 
statistical processing and are presented in the appendices to aid future studies. 
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Figure 7. Convective velocity evaluation. 


D. Concepts 

This section outlines the necessary physical and statistical tools used to check the 
physical model presented in Section III. The model requires a correlation of the undulating 
outer edge of the boundary layer with the region directly above the viscous sublayer. 
Packets of vorticity are formed which migrate from a region called the “production zone.” 
The model predicts that the boundary layer acts as a filter to separate sizes or wave number 
components; that is, the convective velocity becomes a function of wave number. 

Finally,, an asymptotic approach of all sizes of packets to the local mean velocity is 
predicted. Thus, after a finite residence time, the convective velocity will become 
independent of frequency. Methods by which to analyze each of the above hypotheses are 
presented. 
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1. Correlation of the Outer Edge of the Boundary Layer with the Internal Flow . 
The model assumes that the convoluting outer edge of the boundary layer will create strong 
local pressure gradients, which affect what transpires near the wall of the boundary layer. 
However, it is not crucial to know exactly what produces the packet for the major 
contribution of the present report. When the flow near the wall is affected by occurrences in 
the convoluting outer edge, a packet is formed. One may speculate about whether a 
stress-relieving phenomenon [3] may cause localized separation to occur that initiates the 
onset of turbulence. 

A measure of the phenomena proposed is the cross-correlation of the local pressure 
gradients at the two locations in question, such as 


ap \ 

ap \ 

ax,- ) 

' outer 
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i = 1,2; j = 
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inner 


= X-direction; and x 2 = Y-direction. 


Since the local pressure gradients are difficult to measure, one may wish to use local 
accelerations of the fluid particles which are related to the local pressure gradients. Thus, 
the correlation 



outer 



inner 


might be considered. For covariance computations, only the fluctuations need be 
considered : 
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(39) 


Appendix D, equation (D-82), gives the relationship between 
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inner 


31 


and 
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inner . 


The correlation of the velocity fluctuation may be used to evaluate the inner-outer 
correlation described by the model. When only the fluctuating portions of the signals are 
correlated, the product is often referred to as the “covariance.” 

The covariance deals only with the fluctuations. The dc or mean values give no 
insight into convective transport rates, and, thus, only covariance computations will be 
examined for the present applications. “Correlation” will be used interchangeably with 
“covariance,” and, unless otherwise stated, correlation will refer to ac-coupled signals (i.e., 
the means are zero). 

2. Production Zone. The main requirement for this region of the boundary layer is 
that the convective velocity equal the mean velocity; that is, the convective velocity of all 
size structures is the same. 

Evaluation in this zone will be the measure of the convective transport rate for 
various frequencies. If all frequencies travel at the same rate, one can infer the structure’s 
size from the wave number concept. The space time correlation concept discussed in the 
preceding chapter is used to evaluate the convective velocities. 

3. Filtering Action of the Boundary Layer. Filtering of. the boundary layer implies 
an ability of the boundary layer to separate frequencies in some prescribed fashion. The 
physical model predicts that varying frequencies will be convected at different speeds in 
certain portions of the boundary layer. The separation in speed will occur in both the 
vertical and horizontal directions. 

Spatially separated sensors are used to measure the convective speeds of the 
different wave numbers. The convective velocity is used, together with the frequency, to 
determine the wave number. The packet size is related to the wave number. 

4. Convective Velocity Profile Independent of Frequency . The model predicts that 
packets of all sizes will asymptotically approach the mean flow velocity after a finite length 
of time. The example previously given indicates that lengths of the order of 25 to 30 ft (7.6 
to 9.15 m) are required for the convective velocities in a turbulent boundary layer to equal 
the mean velocity. The convective velocities are measured by the spatially separated hot 
wires. 


V. EXPERIMENTAL SETUP AND PROCEDURES 


The experimental portion of this investigation was performed at the Engineering 
Research Center of Colorado State University. Two wind tunnels were used to study the 
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convective velocities: a large meteorological wind tunnel and the entrance portion of a 
small separation tunnel. 

As defined in Section IV, the convective velocity is a time-average quantity, much 
like the concept of mean velocity. Figure 7b is a schematic of the hot wire configuration 
used to define a convective velocity. First, signals from two spatially separated hot wires are 
selectively filtered to obtain desired frequency components. (The wires are separated in the 
X-direction of motion — the mean flow direction.) The signals are then fed into the 
correlator, which produces the signal’s temporal cross-correlation and integrates this value 
for a desired length of time. Integration times of approximately 5 minutes were used for this 
experiment. The cross-correlation is then plotted by an X-Y plotter on graph paper. 

The instrumentaiton involved for this measurement includes the hot wire 
anemometers, filters, a correlator, an X-Y plotter, and a tape transport. 


A. Wind Tunnels 

1. The Meteorological Wind Tunnel . The meteorological wind tunnel test section 
has a cross section that is 6 ft (1.83 m) by 6 ft (1.83 m) and has a 95-ft (29-m) long test 
floor over which the boundary layer may develop. The wind tunnel is of the recirculating 
type with the air speed controlled by a variable speed, variable pitch aircraft propeller. Plate 
and Cermak [68] describe the facility in detail. Figure 1 shows a sketch of the boundary 
layer’s development along the tunnel test floor, and Figure 8 shows a top view of the tunnel. 

Tieleman [66] has shown that damping screens and a 9:1 entrance contraction 
produce extremely low free-stream turbulence, and his results showed that the turbulence 
intensity ranges from 0.00015 at 20 ft/sec (6.1 m/sec) to 0.00027 at 80 ft/sec (24.4 m/sec). 
The boundary layer along the floor of the tunnel is artificially tripped by a saw-toothed 
fence preceded by a 4-ft (1 .22-m) section of gravel (Fig. 7). 

The first 30 ft (9.15 m) of the floor were used for the present experiment. A 
boundary layer thickness of approximately 18 in. (0.457 m) [69] is attained in the 30 ft 
(9.15 m). The measurements were made in the vertical plane along the centerline of the 
tunnel behind the tripped section at distances downstream of 1, 2, 4, 8, 18, and 30 ft 
(0.305, 0.61, 1.22, 2.44, and 9.15 m). The X-axis was taken along the centerline of the 
tunnel floor, the Y-axis is vertical to the floor, and the Z-axis is parallel to the floor and 
perpendicular to the flow direction. 

2. Separation Tunnel. The separation tunnel is of the open return type. A 
honeycomb, which is located in a large circular inlet, is used to straighten the inlet flow and 
break up any large scale disturbances. From the inlet the flow passes through a settling 
chamber and is then accelerated into the test section, where the cross section narrows to 18 
by 18 in. (0.457 m by 0.457 m). The entire test section is made of plexiglass for visual 
observations. A 54-in. (1 .37-m) fan which governs the flow rate by a variable speed motor is 
located downstream of the test section. Figure 9a is a photograph of the flow facility. 
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(b) INSTRUMENTATION 
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5. SCOPE 
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7. FILTERS 

8. 14 CHANNEL TAPE RECORDER 


Figure 9. Wind tunnel and instrumentation. 

The entrance portion of the separation tunnel was used for the initial convective 
velocity measurements. The tests were run at 2 ft, 2 in. (0.66 m) and 4 ft (1 .22 m) from the 
start of the test section. The ceiling of the tunnel was adjusted so as to develop a near zero 
pressure gradient flat-plate boundary layer. The boundary layers were approximately an 
order of magnitude smaller than those of the large wind tunnel. Probes were traversed 
vertically through the boundary layer by an external actuator. The X-axis was parallel to the 
ceiling in the flow direction, the Y-axis was perpendicular to the ceiling in the flow 
direction, and the Z-axis was parallel to the ceiling and perpendicular to the flow direction. 


B. Instrumentation 

1. Hot Wire Probes. The hot wire was used as the basic sensing element for all 
correlation measurements and turbulence measurements. The material for the hot wires was 
80 percent platinum, 20 percent iridium, 0.0004-in. (0.01 -mm) diameter wire. The wires 
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were operated by constant temperature hot wire anemometers designed at Colorado State 
University by Finn and Sandborn [70] . Figure 9b 1 is a photograph of the anemometers. 
Figure 1 0 shows probes of five different designs which were used during the running of the 
experiment. The first probe (Fig. lOal) was used for spatial separation measurements in the 
Z-direction of the inner-outer correlations. One wire was placed close to the wall and the 
other in the intermittent outer edge of the boundary layer. The second probe (Fig. 10a2) 
was a manifold of twelve hot wire elements erected in the Y-Z plane with its major axis in 
the Z-direction. Nine of the twelve wires were in the Z-direction and were used for 
inner-outer correlations as well as spatial correlations. The three remaining wires were in the 
Y-direction for use in space correlations. The next three probes (Fig. 10a) were used in 
traversing the boundary layer for taking mean, turbulent, and convective velocity profiles. 
Each of these probes used two wires: the lead wire being used to make mean velocity and 
turbulence measurements, and the second wire being directly downstream from the lead 
wire but out of the lead wire’s wake. The downstream wire was used to correlate with the 
lead wire to measure a convection time for the flow to traverse from the lead wire to the 
second wire. The length of wires used ranged from 0.05 in. (1.3 mm) to 0.16 in. (4.1 mm). 
The wires ranged from about 3 to 20 ohms cold resistance. Separation distance between the 
leading wire and the lagging wire ranged from 0.002 ft (0.00061 m) to 0.0275 ft (0.00084 
m). Tapered piano wire or jeweler’s broaches were used for the hot wire supports. The piano 
wire or jeweler’s broaches were soldered to conventional copper lead wire, which was 
connected to the anemometer. The hot wire anemometer is comprised of a sensing element 
and a control unit. The output of a constant temperature hot wire anemometer depends on 
the temperature difference between the wire and the local fluid temperature, and the 
velocity of the flowing fluid. There are many different sensing materials which may be used. 

The hot wire anemometer is able to sense a number of fluid properties which may in 
some way transport heat either to or from the sensing element. The temperature of the flow 
was held constant and the pressure fluctuations were very small. It was thus assumed that 
the only first order heat transfer effect was produced by the convection properties of the 
fluid passing over the wire. Therefore, for this case where an overheat ratio of 1:1 is used, 
the “hot” wire will, in general, sense only velocity fluctuations [71], For the convective 
velocity measurements, a calibration of voltage versus velocity is not necessary. Only a 
measure of the ability of the flow to transport its heat transfer qualities is necessary. 

The hot wire is also used for the mean and turbulence measurements. The output 
voltage of the hot wire is a measure of the velocity of the flow. 

To obtain accurate calibration curves, the calibration measurements of the hot wires 
were taken in the freerstream section of the wind tunnel against a calibrated pitot tube. An 
averaging circuit, included as an integral part of the anemometer housing, was used to obtain 
an integrated averaged voltage. 

Each hot wire was annealed before its calibration and use. That is, the wire was 
heated until red hot, and was allowed to “cook” until the hot wire characteristics became 
stabilized. The hot wires were calibrated before and after each run. A typical calibration 
curve is given in Figure 1 1. The slope of the calibration curve 3E/9U is the sensitivity of the 
hot wire to velocity. Velocity rms values were calculated by measuring the rms and the 
mean value of the signal, with the local sensitivity, dE/dU , being determined graphically at 


36 



(a) VARIOUS HOT WIRES 

1. ORIGINAL INNER -OUTER PROBE 

2. 12 WIRE MANIFOLD FOR INNER- 
OUTER CORRELATIONS AND 
SPACE CORRELATIONS 

3. 2 WIRE PROBE FOR CONVECTIVE 
VELOCITIES 

4. 2 WIRE PROBE FOR CONVECTIVE 
VELOCITIES 

5. 2 WIRE PROBE FOR CONVECTIVE 
VELOCITIES 


( b) PRINCETON CORRELATOR 

1. GAIN CONTROLS 

2. COUPLING CONTROL 

3. SIGNAL DELAY CONTROL 

4. READOUT RATE 

5. INPUT CHANNELS 

6. OUTPUT OF TIME BASE 

7. OUTPUT OF CORRELATION 




Figure 10. Instrumentation. 

the mean value of the signal, E. Multiplying 8U/3E with the rms of voltage, >/?”, gives the 

rms of the velocity, . An error could arise in this evaluation if the turbulence intensity 
becomes too great, although this is not a problem for turbulence intensities below 
approximately 15 to 20 percent [71]. High turbulent intensities are encountered near the 
boundary, and these measurements are therefore less accurate. The sensitivity of a hot wire 
to flow direction and an empirical solution for the hot wire evaluation are given by Zoric 
[67]. Sources of error such as solid boundary effects, wire length, velocity gradient, and 
turbulent intensity are covered by Tieleman [66] and Zoric [67] . 

For convective velocity measurements, a new source of error may be a phase or 
time-delay mismatch of the anemometer electronics. The convective velocity measurement 
relies on the fact that time-delay separation of leading and lagging wires is dependent only 
upon the convective flow velocity. The two hot wire anemometers used in the research were 
checked for differences in time delay or phase shift. By placing both of the hot wire sensing 
elements very close together in a turbulent flow, and cross-correlating the output of the two 
systems, any time-delay errors would be seen as a difference between the auto- and 
cross-correlations. Figure 12 shows the result of one such test. The units were found to be 
identical for the purposes of this experiment. 
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2. Pitot-Static Tube . Mean velocity was measured by a standard 1/8-in. (3.18-mm) 
diameter pitot-static tube. The pitot tube used for mean velocity measurements was 
calibrated by comparing it with a laboratory standard pitot tube. The' standard had 
previously been calibrated by placing it on a whirling arm, the length of which is known. 
With this information, along with the number of revolutions made per unit time and 
correcting for swirl effects, the laboratory standard was calibrated. The pitot pressure from 
the laboratory standard was read on a water, U-tube micromanometer, as was the case for 
the test pitot tube. The micromanometer has been calibrated from laboratory standard 
micromanometers. The possible sources of error for the mean velocity measurement are 
covered by Zoric [67] . They were found to be negligible for this experiment. 

3. Actuator. Two types of actuators were used for precision placing of the sensors 
in the vertical direction of the boundary layer. The first was a stand-type model, which was 
set directly on the floor of the meteorological wind tunnel. A small electric motor with a 
screw-type gear was used to change the probe elevation. A variable potentiometer which 
changed resistance as the carriage moved up or down was used to determine the elevation of 
the probes. The second type of actuator extended a rod through the ceiling of the 
separation wind tunnel. Variations in the external length of the rod gave a direct measure of 
the probe location in the flow. This actuator was also controlled with an electric motor and 
a resistance potentiometer with an emf output to measure the elevation of the probe. The 
first actuator, used in the meteorological tunnel, had a range of approximately 18 in (0.46 
m) while the second actuator had a range of about 6 in. (0.15 m). In both cases initial 
conditions of probe placement were determined with a standard length scale, within an 
accuracy of 0.005 in. (0.13 mm). 

4. Electronic Filters. The filters used in this study were continuously adjustable 
from 0.2 Hz to 20 kHz in either a low-pass or high-pass mode. Filter attenuation and phase 
shift as a function of frequency were found to be identical for the two filter banks 
employed. Attenuation of the filter units shown in Figure 9b, was 24 dB per octave. 

Since filters may cause nonlinear phase shifts as a function of frequency, two 
electronic filtering systems that are not identical would produce what appears to be a time 
delay between two identical signals after passing through the filters. Thus, a calibration and 
check procedure was used to determine phase differences as a function of frequency for the 
filters used in this experiment. All possible types of auto- and cross-correlations of filtered 
hot wire signals were used to check for phase differences. The same signal was fed through 
each filter and then correlated. Similarly, identical sine waves were used. Both cases showed 
the filters to be identical as far as the present study was concerned. 

Figures 13 and 14 show the phase checks using a hot wire signal as the emf. 

5. Time Correlator . For auto- and cross-correlation of hot wire signals, the 
Princeton Applied Research Correlation function computer model 101 was used (Fig. 10). 
The correlator is a hybrid computer which uses both analog and digital techniques. The 
correlator shifts one signal in time, multiplies the two signals together, and integrates. It 
performs the operation for 100 discrete points each separated by a time shift of Ar . The 
time shift At between successive points may be selected. The correlation may be examined 
from 0 to 100 At in steps of At . A precomputation period may also be selected, so that 
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Figure 13. Phase check of filters. 
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Figure 14. Correlation phase check of filters using typical hot wire signal. 


the initial time delay is not 0, but j( 1 00 At), where j may be selected from 1 to 10. The 
term 100 A r may be selected from 0.0001 to 10 sec. Readout circuitry allows the 
continuous monitoring of the function during operation. The. correlation function may be 
either viewed on an oscilloscope or plotted on an X-Y plotter or strip chart recorder. The 
accuracy of the unit is within 1 percent of the ideal function. Both ac and dc couplings are 
available. 

The correlator was checked for phase shift and time delay within its own dual 
circuitry by cross-correlating a sine wave from a calibrated function generator. This check is 
good for both phase shift and time delay calibration. The magnitude of. the correlator 
compared with a calibrated rms meter was found to be in calibration. 

6. Tape Transport. The tape transport used in this experiment was a 14-channel 
standard I-rig (Fig. 9b) recorder. The recorded signals were FM modulated. The carrier 
frequencies were calibrated with a special commercial unit. The 14 channels are split into 
two sets of heads, one set for odd numbered channels and one set for even numbered 
channels. Only odd channels, or only even channels, were used when recording a run to 
avoid time-displacement errors caused by the physical separation of the heads. Interchannel 
displacement errors were checked by recording various frequency sine waves and correlating 
the outputs. Any off-zero peak of the correlation would indicate a time-displacement error 
in the recorder. For frequencies less than 100 kHz, no error could be detected. Frequencies 
of interest for this experiment were below 2 kHz. Dynamic interchannel time placement 
errors in the recorder were also checked and found to be inconsequential. 

7. Other Instrumentation. The ordinary types of instrumentation which are 
common will not be described in detail, but rather are listed so that the reader will have a 
feeling for the type of instrumentation needed for this type of experiment. They are 
oscilloscopes, function generators, rms meters, digital voltmeters, and X-Y plotters. The 
type of transport described above is also desirable, so that the identical signal may be 
analyzed for all computations, some of which must be done recursively. 

The time base of the correlator was used to calibrate the X-component of the X-Y 
plotter, and a standard cell was used for amplitude calibration. 


C. Convective Velocity Computations 

The convective velocity was determined by the following procedure: 

1 . Orient two spatially separated hot wires in the direction of the mean flow. 

2. Measure the magnitude of the hot wires’ separation. 

3. Cross-correlate the output from the two spatially separated sensors. 

4. DetermineThe time delay for which the cross-correlation is a maximum. 

5. Now, dividing the separation distance by the time delay for which the correlation 
was a maximum gives the convective velocity. We see that simultaneous filtering of phase 
matched systems may be used to describe convective velocities of various frequencies. 
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D. Inner-Outer Correlation Measurements 

The velocity fluctuations in the intermittent outer region were correlated with 
velocity fluctuations near the wall using two hot wire sensing elements in the respective 
locations. Signals were put directly into the correlator, which was described previously. The 
normalized correlation is used to give a relative measure of the similarity of the two signals, 
and it gives some insight into a possible mechanism for the triggering of turbulent 
production near the wall. A time delay of the correlation other than zero would indicate the 
length of time for the production process to acknowledge a stimulus produced at the outer 
edge. 


I r 

E. Data Reduction of Time Series 

Because of the length of this subject it has been placed in Appendix D. General 
signal analysis is covered, as well as some extensions which give insight into methods of 
dealing with convective velocities in Fourier transformation planes. Restrictions imposed 
upon correlations by either finite time series lengths or certain degree's of nonstationarity of 
the time series are examined. The effect of space and frequency weighting functions on the 
normalized cross-covariance is analyzed, and finally the space-time covariance is viewed as to 
its relationship with the probability density distribution of the original, time series. 


VI. EXPERIMENTAL RESULTS AND DISCUSSION 


The theoretical model presented in Section III predicts the occurrence of several 
phenomena, some of which have been confirmed by previous investigators. 

The objective of this section is to present experimental evidence confirming the 
physical model presented earlier. Experimental evidence shows correlations of the 
intermittent outer edge with flow near the wall. Data confirming the existence of a 
formation zone where the convective velocity is independent of wave number are presented. 
The convective velocity is seen to be a function of wave number throughout most of the 
boundary layer. The convective velocity is shown to approach the mean velocity in the 
outer portion of the boundary layer at a distance downstream, as predicted by the model. 
Several special cases concerning the convective velocities will also be examined. 

A. Outer-Inner Correlation 

The correlation between the intermittent outer edge of the boundary layer with 
flow near the wall was evaluated in the small wind tunnel previously described. The 
boundary layer thickness was about 3 in. (0.076 m) with a free-stream velocity of 
approximately 30 ft/sec (9.15 m/sec). The basic principle involved was to determine if a 
correlation existed between the velocity fluctuations of the undulating outer edge of the 
boundary layer and the turbulence near the wall. 
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The outer-inner correlation was performed by placing a hot wire in the highly 
intermittent outer region of the boundary layer and one very near the surface. Both wires 
were oriented normal to the flow and parallel to the surface. The size and type of wires used 
are described in paragraph V.B.l. A sketch of the setup is given in Figure 15b. The signals 
were ac-coupled and cross-correlated. The normalized cross-covariance, was then computed. 
The experimental values of the correlation obtained were about 0.1. Figure 15a shows a 
typical experimental.' curve. It was always found that a time delay existed between the two 
sensors for the peak of the correlation; that is, the maximum correlation occurred when the 
signal from the wire near the wall was' delayed slightly in time before correlating it with the 
wire in the intermittent outer region. In an attempt to see if it was possible to get a 
maximum of the correlation to occur nearer zero time delay, the lower sensor (sensor y in 
Fig. 15b) was slowly moved upstream. Moving the probe upstream so that the lower probe 
was either directly under or upstream of the outer probe resulted in the complete 
elimination of a measurable correlation. A measurable correlation could only be obtained 
when the sensor near the surface was slightly downstream of the outer probe. This result 
seems to give support to the assumption that disturbances are transmitted from the outside 
of the boundary layer to the inside, and not vice versa. The vanishing of the correlation, 
when the inner probe is moved upstream of the outer probe, gives some confidence that the 
correlation did not result from large scale circular structures that pass through both probes 
simultaneously. 

Other configurations of wire geometry were also used for this correlation; for 
example, the outer wire parallel to the mean flow and the inner wire normal to both the 
mean flow and floor. The configuration where both wires were perpendicular to the flow 
produced the same result as that above. The manifold-type of probe (Fig. 10b) was also used 
with the same results. These measurements are in agreement with the physical model 
presented at the beginning of this, report. Laufer and Narayanan [59] have shown that 
processes occurring near the viscous sublayer scaled with outer flow parameters. . 

The experimental test run, correlating the inner and outer edge of the boundary 
layer, tempts one to speculate about the role of the outer edge in the onset of turbulence 
near the wall. It may be inferred that information is fed into the near surface region of the 
boundary layer from the undulating intermittent outer edge. 


B. The Formation Zone 

• • ■ 1 

The physical model required the assumption that there is a region near the surface 
where the convective velocity is identically equal to that of the local mean velocity. 
Convective and mean velocity measurements were made simultaneously across the boundary 
layefr to determine if such a condition did exist. 

( The convective velocity was calculated using two spatially separated hot wires. The 
wires were spatially separated in the direction of the mean motion. The lead hot wire was 
used to measure mean velocity as well as the turbulent velocities. The cross-correlation of 
the ac-coupled lead wire with the ac-coupled second wire was evaluated. A maximum 
correlation at some finite time delay was obtained. (Figure 16 is a typical example of the 
correlations that were found.) The spatial separation of the wires divided by the time delay 
at which the maximum correlation occurred was then considered a measure of the 
convective velocity; that is, 
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Figure 1 6. Correlation curves with no filtering in the outer region. 
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Figure 6 gives a graphical picture of this process, where C y is the convective velocity in the 
direction of the spatial separation, AW , and r m is the time delay at which the 

cross-correlation is a maximum. Section IV describes the convective velocity "concept in 
more detail. 

Figures 17 and 18 show the experimental results obtained across two different 
boundary layers. (Note: Two different wire separations were used for these two sets of 
measurements.) Both of these flat plate velocity profiles give a region where the mean and 
convective velocities are equivalent. It appears that this may be more likely a region than a 
point, where the two profdes are very nearly equal. The location of the region corresponds 
closely with that observed by Kline et al. [2, 33] in the visualization studies. The convective 
velocity is less than the mean velocity in the outer portion of the boundary layer. At 4 the 
same time, the convective velocity is higher than the mean velocity profile below the 
formation zone. 

The problems associated with one wire following downstream of another must be 
evaluated. One must be sure that the trailing wire is not affected by the wake of the 
upstream wire or probe. This was checked by moving the front probe slightly, and seeing if 
there was an effect on the dc level of the second wire. That is, if the second wire was in a 
wake flow caused by the leading wire or probe, the mean velocity recorded by the second 
wire would be less than if the front wire or probe was not there. The fact that the 
convective velocity exceeds the mean velocity indicates that the convective velocity is not a 
probe-interference phenomenon and gives some assurance that the second probe is not in 
the wake of the first. Closely spaced, large wires would be more prone to interference errors. 
The rear wire was also shifted up and down to see if a velocity defect region existed. There 
was no indication that the second wire was in the wake of the first. 

i 

The cross-correlations of the longitudinal velocity fluctuations across the boundary 
layer are seen in Figures 16, 19, and 20. It is interesting to note the sudden decrease in 
amplitude of the correlation near the wall. The decrease is in part due to the fact that the 
high frequency spectral components are dissipated in the viscous subregion. 

Appendix D, paragraph F, shows the relationship between the probability density 
distribution and the covariance computation. Equation (D-93) shows that, since the 
covariance is a squared function of the perturbation size and a linear function of time, the 
larger excursions from the mean are the most heavily weighted. It is also found that if the 
turbulent velocities continue to propagate with the velocity indicated at the first sensor, 
equal perturbations above and below the mean velocity will create correlations of equal 
magnitude at time delays equally spaced about the time delay that indicated the mean 
velocity. ■ . >' • - : • ' ' ■ r . 

It is anticipated that, a symmetrical probability density distribution will produce a 
correlation peak that will indicate the mean velocity. Figure 21 shows the probability 
density distribution and probability distribution at specific points in the flat plate boundary 
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Figure 17. Mean and convective velocity versus elevation, AW = 0.0233 ft (7.1 mm). 
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Figure 18. Mean and convective velocity versus elevation, AW = 0.0183 ft (5-.58 mm). 


CORRELATION CORRELATION 



At 



Figure 19. Correlation curves with no filtering. 

layer, which was about 2 in. (5.08 cm) thick at the point of measurement. Figures 21, a and 
b, show the distributions near the outer edge of the boundary layer. The obvious skewing 
toward the negative side indicates that the convective velocities will be substantially lower 
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At 

Figure 20. Correlation curves with no filtering near the surface. 

than the mean velocity at these locations. This indication of a lower convective velocity 
agrees with the actually measured convective velocities. Figures 21, c through f, show the 
probability densities closer to the wall. The distributions are skewed over a large part of the 
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PROBABILITY DENSITY DISTRIBUTION 

— PROBABILITY DISTRIBUTION 



- e + - e + 


(a) 1.900 in. (4.83 cm) FROM PLATE (b) 1.100 in. (2.79 cm) FROM PLATE 



- e + - e + 


(e) 0.100 in. (0.254 cm) FROM PLATE (f) 0.025 in. (0.0635 cm) FROM PLATE 

Figure 21. Probability distributions in the flat plate boundary layer, 26 in. (0.66 m) 
from start of test section, free stream 40 ft/sec (12.2 m/sec). 

boundary layer. Only very near the wall does the probability distribution appear 
symmetrical. A symmetrical probability distribution would indicate that the convective 
velocity is almost equal to the mean velocity (Appendix D, paragraph F). 
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As noted in Appendix D, paragraph F, if the probability distribution indicates a 
skewing toward the positive or negative velocity fluctuations, the convective velocity will 
likewise be correspondingly higher or lower than the mean velocity. The correlation 
function is built up by all the velocity contributions of the probability distribution. The 
wings of the probability distribution are the most heavily weighted. A convective velocity 
less than the mean would give a time delay greater than the time delay for the mean 
velocity; a convective velocity greater than the mean would give a time delay less than that 
for the mean velocity. Appendix D, paragraph F, gives the contribution produced by each 
turbulent velocity acting over some Finite time interval. Thus, a skewing of the probability 
distribution is propagated into a skewing of the shape of the correlation function. Figures 
16, 19, and 20 are photographs of the correlation function for various Y-distances. Note 
that the correlation curve skews to the right at the outer edge of the boundary layer and 
becomes more symmetrical nearer the surface. The shapes of the probability distributions 
and correlation functions are consistent with the proposed physical model. The symmetries 
of the correlation function and probability density distribution are consistent with the 
production zone of the model. Here the convective velocity is equal to the mean velocity, so 
that symmetrical correlations and probability density distributions are observed. 


C. Convective Velocity Profiles as Functions of Frequency 

This section presents experimental data to support the physical model’s hypothesis 
that the boundary layer will frequency-filter the packets of vorticity that migrate from the 
production zone; that is, the packets will be separated according to frequency. Also, some 
support is given to justify associating frequency with size. 

The analysis in Section IV required that the boundary layer local velocity separate 
the various scales as they migrate from the production zone. In general, only at the point of 
origin, where all scales are moving with the mean velocity, should the convective velocity be 
independent of frequency. The special case when the convective velocities have reached 
equilibrium with the mean flow will be discussed later. 

Convective velocity measurements were taken across a flat plate, zero pressure 
gradient boundary layer using various wire separations. Willmarth et al. [51, 52] noticed 
that the convective velocity of pressure fluctuations on the floor of a flat plate boundary 
layer were functions of sensor separation. The measurements taken during this experiment 
also showed that for very closely spaced sensors, the convective velocities are a function of 
sensor spacing. The convective velocity for closely spaced sensors is lower than convective 
velocities measured with larger separations for the same frequencies. 

Figure 22 shows the measured mean velocity profiles along with convective velocity 
profiles for various frequencies. Phase-matched filters were used together with the two hot 
wires, as indicated in paragraph VI. A, to obtain Figure 22. The normal dropoff of the 
turbulent spectrum acted as a low pass filter, while electronic filtering was used for the high 
pass filtering- This combination of electronic filtering and turbulent spectral rolloff acts as a 
notch filter, which allowed the evaluation of convective velocity for narrow bandwidths. 
Figure 22 represents convective velocities measured with different separations between the 
sensors. In all cases the farther from the production the measurement is made, the greater 
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0.15 NUMBERS INDICATE HIGH 
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Figure 22. Mean and convection velocity versus elevation for various AW. 

the spread is of the convective velocities, and, as the profiles get nearer the surface, the 
spread of the convective velocities decreases and begins to collapse to a single velocity. 
Another way to view this problem is to plot the ratio of the convective velocity at a given 
frequency to the convective velocity with notch filtering versus wave number. The wave 
number is computed by the formula, 


K f = 



(41) 
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where is the convective velocity at frequency f (V c will sometimes be expressed 
as Vj- for brevity). Figure 23 shows the variation of the data of Figure 22 plotted in 


convective velocity versus wave number form. In every case the profile of the convective 
velocity versus wave number which is nearly one for all wave numbers 1 is that of the 
measurement nearest the production zone (Tables 1 through 4). The spreading or frequency 
filtering of the convective velocities by the boundary layer was predicted by the physical 
model. The spreading of all frequencies from a single point supports the belief that this is a 
formation zone. 


Figure 24 shows. measurements in the region both above and below the formation 
zone in a very large boundary layer. The large boundary layer showed the same effects as 
the small boundary layer. As before, spreading occurs whenever the. measurement is made 
away from the production zone (Table 5); that is, above the production zone the higher 
frequencies moved faster than the lower frequencies, and below the production zone the 
higher frequencies move sjower than the lower frequencies. However, in both cases the 
higher frequencies more closely approach the mean velocity than do the lower frequencies. 
In general, the convective velocities above the production zone are slower than the' mean 
velocity; the convective velocities, below the production zone are higher than the mean 
velocity. The mean velocity obviously is the driving or forcing function on any packet or 
structure in the flow. The relationship between convective and mean velocities, as indicated 
above, is consistent with the physical model presented. Figure 25 shows a typical set of 
cross-correlations for convective velocities in the outer portion of a boundary layer. Note 
that the skewing is consistent with the analysis of Appendix D, paragraph F. The 
progression of higher frequencies toward the mean velocity time delay is consistent with the 
physical model. • , 

The coalescing of all convective velocities at a particular point makes it difficult to 
accept models which require large eddies with characteristically high velocities protruding 
into the lower layers of the boundary layer. 


Figures 22 and 23 and Tables 1 through 4 represent four separate evaluations of 
convective velocity with various wire separations. Figures 26 and 27 present the above data 

■ ■. aw 

in the .form of V c ^/V m versus K f^T • Although V c , V m , Kf , and AW have 

previously been defined, has not in itself been examined. If one considers a time 

function with an angular frequency of 2nf convecting with a' velocity .V c ~. ? the number of 

cycles that the function performs while traversing a distance AW is 



Vc.o/ V c ,f 



Figure 23. Convective velocity versus wave number. 

This dimensionless number is similar to a dimensionless frequency or Strouhal number. 
However, here the parameter represents the number of oscillations a particular frequency 
will make between the first and second sensor. Figure 26 shows the individual test results. 
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TABLE 1. TEST I CONVECTIVE VELOCITY PARAMETERS 


Run No. 1 

V m = 40.4 ft/sec (12.3 m/sec) ' 

AW = 0.01 167 ft. (3.56 mm). ... Elevation = 1.968 in. (4.99 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

Vf/Vo 

. K f 

Ml 


All 

27.8 

.1.00 

0.00 

n 

0.0000 

1' 

27.8 

1.00 

, 0.23 

■HU 

0.0004 

10 

27.8 

1.00 

2.26 

0.688 

0.0042 

100 

27.8 

.1.00 ' 

22.60 

0.688 

0.0419 

200 

27.8 

1.00 

45.20 

0.688 

0.0840 

300 

27.8 

1.00 

'67.80 

0.688 

0.1259 

400 

' 28.6 

1.029 

87.88 

0.708 

0.1632 

500 

29.15 

1.049 

107.8 

0.722 

0.2002 

600 

30.4 

1.094 

124.0 

0.752 

0.2303 

700 

30.7 

1.104 

143.3 

0.760 

0.2261 

800 

31.35 

1.128 

160.3 

0.776 

0.2978 


Run No. 2 ...... 

V m = 39.6 ft/sec (12.1 m/sec) Elevation = 1 .788 in. (4.55 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

V f /Vo ' 

K f 

V V m 

K f AW/27r, 

cycles 

All 

27.8 

1.000 

0.000 

0.702 

0.0000 

1 

' 27.8 

1.000 

0.226 

0.702 

0.0004 

10 

27.8 

1.000 

: 2.260 

0.702 

0.0042 

100 

27.8 

1.000 

22.60 

0.702 

0.0420 

200 

27.8 

1.000 

45.20 

0.702 

0.0840 

300 

27.8 

1.000 

67.80 

0.702 

0.1260 

400 

• 27.8 

1.000 

90.40 

0.702 

0.1680 

500 

• 28.6 

1.029 

109.8 

0.722 

0.2040 

600 

29.2 

1.050 

129.1 

0.737 

0.2398 

" 700 

30.1 

1.083 

146.1 

0.760 

0.2714 

800 

30.7 

1.104 

163.7 

0.775 

0.3041 

900 

31.7 

1.140 

178.4 

0.801 

0.3313 

1000 

31.7 

1.140 

198.2 

0.801 

0.3681 
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TABLE 1. (Continued) 


Run No. 3 

V m = 35.9 ft/sec (1 1.0 m/sec) Elevation = 1.053 in. (2.67 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

"t, 

> 

K f 

V V m 

K f AW/27r, 

cycles 

All 

24.3 


0.000 

0.67.7 


1 

24.3 


0.260 

0.677 


10 

24.3 

. 

2.580 

0.677 


100 

24.3 


25.86 

0.677 

0.0480 

200 

24.3 


51.71 

0.677 

0.0960 

300 

25.2 

■V 

74.80 

0.702 

0.1389 

400 

25.3 

1.041 

99.34 

0.705 

0.1845 

500 

25.8 

1.062 

121.7 

0.719 . 

0.2262 

600 

26.3 

1.082 

143.3 

mEam 

0.2662 

700 

27.0 

1.111 

162.9 


0.3026 


Run No. 4 

V m = 26.4 ft/sec (8.05 m/sec) Elevation = 0.368 in. (0.935 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

> 

Kf 

V f/ V m 

K f AW/2ir, 

cycles 

All 


HU 

0.000 

0.773 

0.0000 

1 


mSSm 

0.308 

0.773 

0.0006 

10 

20.4 

Bees 

3.080 

0.773 

0.0057 

100 

20.4 

mm 

30.80 

0.773 

0.0572 

1/ 


■Ea 

61.60 

0.773 

0.1144 

1 

20.6 


91.50 

0.780 

0.1699 

400 

20.8 


120.8 

0.788 

0.2244 

500 

21.0 

mm 

149.6 

0.795 

0.2779 

600 

21.5 

mam 

175.3 

0.814 

0.3257 

700 

21.6 

1.059 

203.6 

0.818 

0.3782 

800 

21.8 

1.069 

230.6 

0.825 

0.4283 
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TABLE 1 . (Concluded) 


Run No. 5 

V m = 18.7 ft/sec (5.7 m/sec) Elevation = 0.050 in. (0.127 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

> 

*f 

v f/v m 

K f A W/27:, 
cycles 

All 

17.2 


0.00 

0.920 

0.0Q00 

1 

. 17.2 


0.37 

0.920 

0.0007 

10 

17.2 


3.65 

0.920 

0.0068 

100 

17.2 

■ 

36.53 

0.920 

0.0678 

200 

17.2 


73.06 

0.920 

0.1357 

300 

17.2 

1.000 

109.6 

0.920 

0.2035 

400 

• 17.2 

1.000 

145.1 

0.920 

0.2695 

500 

17.2 

1.000 

172.6 

0.920 

0.3207 

600 

17.25 

1.003 

218.6 

0.922 

0.4059 

700 

17.35 

1.009 

253.5 

0.928 

0.4708 

800 

-17.35 

1.009 

289.7 

0.928 

0.5381 


TABLE 2. TEST II CONVECTIVE VELOCITY PARAMETERS 


Run No. 1 

AW = 0.0271 ft (8.25 mm) 

V m = 41.0 ft/sec (12.5 m/sec) Elevation = 3.670 in. (9.33 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

V f /Vo 

Kf 

v f/ v ra 

K f AW/27r, 

cycles 

All 

27.4 

Bi 

0.00 

0.668 

0.0000 

10 

27.4 : 


2.29 

0.668 

0.0099 

100 

! 27.4 

HESS 

22:93 

0.668 

0.0989 

200 

• 27.4 

MEM 

45.86 

0.668 

0.1978 

300 

28.0 

iUQSt 

67.32 

0.683 

0.2904 

400 

28.6 

1.04 

87.87 

0.698 

0.3790 

■ 500 

29.4 

1.07 

106.8 

0.717 

0.4609 

: 600 

S 30.4 

1.11 

124.0 

0.741 

0.5349 

700 

31.0 

1.13 . 

141.9 

0.756 

0.6119 

800 

31.5 

1.15 

159.6 

0.768 

0.6882 

900 

31.9 

1.16 

177.3 

0.778 

0.7646 

1000 

31.9 

1.16 

196.9 

0.778 

0.8496 
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TABLE 2. (Continued) 


Run No. 2 

V m = 40.7 ft/sec (12.4 m/sec) Elevation = 2.62 in. (6.65 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X"3.28) 

< 

< 

o 

Kf 

Vf/V m 

K f AW/2rr, 

cycles 

All 

26.0 

mm 

0.00 

mm 

0.0000 

10 

26.0 

ilRB 

2.24 

— 

0.0104 

100 

26.0 

Mat 

24.17 

0.639 

0.1042 

200' • 

27.1 

1.04 

46.37 

0.639 

0.1200 

300 

27.6 

1.06 

68.30 

0.666 

0.2946 

400 

28. 2 : 

1.08 

89.12 

0.693 

0.3844 

500 

29.0 

1.12 

108.3 

0.713 

0.4672 

600 

30.0 

1.15 

125.7 

0.737 

0.5420 

700 

30.8 

1.18 

142.8 

0.757 

0.6159 

800 

32.2 

1.24 

156.1 

0.791 

0.6733 

900 

32.5 

1.25 

174.0 

0.799 

0.7505 

1000 

33.0 

1.27 

190.4 

0.811 

0.8212 


Run No. 3 

V m = 40.5 ft/sec (12.4 m/sec) Elevation = 1.820 in. (4.63 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

O 

> 

> 

Kf 

v f/ v m 

K f AW/27r, 

cycles 

All 

25.8 

mm 

0.00 

mam 


100 

26.1 

mm 

24.07 

mBm 


200 

26.3 

wEm 

47.78 



: 300 

27.1 


69.56 

Wmm 

0.3000 

400 

27.6 

1.07 

91.06 

m 

0.3928 

500 

28.2 

1.09 

111.4 

0.696 

0.4805 

600 

29.4! 

1.14 

128.2 

0.726 

0.5531 

700 

30.0 

1.16 

146.6 

0.741 

0.6323 

: ’800 

30.8 

1.19 

163.2 

0.760 

0.7039 

900 

31.2 

1.21 

181.2 

0.770 

0.7817 

1000 

31.6 

1.22 

198.8 

0.780 

0.8576 
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TABLE 2. (Continued) 


Run No. 4 

V m = 36.5 ft/sec (1 L2 m/sec) Elevation = 1.220 in. (3.1 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

> 

Kf 

>f/ V m 


All 

24.2 

mm 

0.OO 

mm 


10 

24.3 

■MS 

2.59 


0.0111 


24.8 

IHjEfl 

25.34 

0.666 

0.1117 


25.0 

BES 

50.27 

0.679 

0.2168 

I V 

25.8 

1.07 

73.06 

0.685 

0.3151 

400 

26.2 

1.08 

95.93 

0.707 

0:4138 

500 

27.4 

1.13 

114.7 

0.718 

0.4945 

600 

27.6 

1.14 

136.6 

0.751 

0.5891 

700 

28.2 

1.17 

155.9 

0.756 

0.6727 

800 

28.4 

1.17 

176.9 

0.773 

0.7634 

900 

28.6 

1.18 

197.7 

0.778 

0.8528 

1000 

29.0 

1.20 

216.7 

0.795 

0.9345 


Run No. 5 . 

V m = 30.85 ft/sec (9.4 m/sec) Elevation = 0.635 in. (1 .61 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

v f /v 0 

. K f . 

v f/ v m 


All 

21.7 

1.00 

0.00 


■n 

10 

21.7 

1.00 

2.90 



100 

21.7 

1.00 

28.95 

1 

ins 

200 

22.0 

1.01 

57.12 

0.713 


300 

22.5 

1.04 

83.78 

0.729 

0.362 

400 

23.6 

1.09 

106.5 

p.765 

0.457 , 

500 

23.8 

1.10 

132.0 

0.771 

0.570 

600 

24.2 

1.11 

155.78 

0.784 

0.671 

700 

24.4 

1.12 

180.3 

0.7,91 

0.787 

800 

24.4 

1.12 

206.0 


0.890 

900 

24.5 

1.13 

230.8 

0.794 

0.994 

1000 

24.5 

1.13 

» 

256.5 

,0.794 

1.106 
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TABLE 2. (Concluded) 


Run No. 6 

V m = 25.5 ft/sec (7.8 m/sec) Elevation = 0.240 in. (0.61 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

V f /Vo 

K f 

v f/ v m 


All 

17.1 

1.00 

0.00 

0.671 


10 

17.1 

1.00 

3.67 

0.671 

0.016 

100 

17.1 - 

1.00 

36.74 

0.671 

0.158 

200 

17.1 

1.00 

73.49 

0.671 

0.317 

3 00 

17.1 

1.00 

110.0 

0.671 

0.475 

400 

17.2 

1.01 

146.1 

0.671 

0.631 

500 

17.2 

1.01 

182.7 

0.671 

0.788 

600 

17.3 

1.01 

212.1 

0.671 

0.915 

700 

17.3 

1.01 

254.2 

0.671 

1.098 

800 

17.3 

1.01 

290.6 

0.671 

1.253 


TABLE 3. TEST III CONVECTIVE VELOCITY PARAMETERS 


Run No. 1 

AW = 0.002083 ft. (0.635 mm) 

V m = 35.7 ft/sec (10.9 m/sec) Elevation = 1 .105 in. (2.81 cm) 


Convective 

Frequency 

' Convective 
Velocity, 
fps 

(m/sec X 3.28) 

< 

< 

O 

% 

V f/ V m 


AH 

7.44 

1.000 

0.00 

SE9 

0.0000 

10 

7.44 

1.000 

8.45 

WmM 

0.0028 

100 

7.44 

1.000 

84.45 

WBm 

0.0280 

200 

8.0.4 

1.081 

156.3 


0.0530 

300 

8.67 

1.165 

217.4 

mtm 

0.0720 

400 

10.4 

1.398 

241.7 

0.291 

0.0800 

500 

13.0 

1.747 

241.7 

0.364 

0.0800 
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TABLE 3. (Continued) 


Run No. 2 

V m = 30.5 ft/sec (9.3 m/sec) Elevation = 0.525 ini (1.33 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/ sec X 3.28) 

'V f /V 0 

K f 

■ V f/ V m 


All 

■ 8.0 


0.000 

0.262 


10 

8.0 

mm 

7.850 

0.262 


100 

8.0 

■ms 

.78.54 

0.262 

0.0260 

200 

5 • 8.0 

■fri 

157.1 

0.262 

0.0520 

300 

8.5 


221.8 

0.278 

0.0732 

400 

9.9 

■tfl 

253.9 

0.324 

0.0838 

500 

11.6 

WSm 

270.8 

0.380 

0.0895 


Run No. 3 

V m = 26.5 ft/sec (8.1 m/sec) Elevation = 0.265 in (0.673 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

> 

Kf 

v f/v m 


All 

7.4 

1.000 

0.000 

0.279 

0.0000 

10 

7.4 

1.000 

8.490 

0.279 

0.0028 


7.4- 

1.000 

85.91 

0.279 - 

0.0281 


7.4 

1.000 

169.8 

0.279 

0.0562 


8:3 

1.122 

227.1 

0.313 

0.0751 


9.9 

1.338 

253.9 

0.374 . 

0.0838 


*• 11.6 

1.568 

270.8 

0.437 

0.0895 


Run No. 4 . 

V m = 21.2 ft/sec (6.47 m/sec) Elevation = 0.080 in. (0.203 cm) 


Convective 

Frequency 

Convective 
Velocity, .. 

. fps 

(rn/sec X 3.28) . 

.. Vf/Vo . 

-Kf 

■ 


All 

5.3 

1.000 

0.00 

■ms 

0.0000 

10 

5.3 

1.000 

11.86 


0.0039 

100 

5.3 

1.000 

118.55 

BH 

0.0393 

200 

5.3 

1.000 

237.1 


0.0785 

300 

6.2 

1.170 

303.7 

wmtm 

0.1005 

400 

6.9 

1.302 

364.2 

mm 

0.1207 
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TABLE 3. (Concluded) 


Run No. 5 1 

V m = 15.7 ft/sec (4.8 m/sec) Elevation = 0.020 in. (0.0508 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

> 

K f 

V f/ V m 


All 

5.5 

mm 



0.0000 

! 10 

5.5 

BK'tw 

11.42 


0.0038 

100 

5.5 

mm 

114.2 

0.350 

0.0378 


5.5 

■Bl 

228.5 ■ 

0.350 : 

0.0756 


5.8 

■HI 

324.9 

0.370 

0.1077 


TABLE 4. TEST IV CONVECTIVE VELOCITY PARAMETERS 


Run No. 1 

AW = 0.0050 ft (15.3 mm) 

V m = 41.2 ft/sec (1 2.6 m/sec) Elevation = 2.708 in. (6.9 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

Vf/Vo 

K f 

V f/ V m 


All 

2.08 

1.000 

0.00 

0.051 

0.0000 

.10 

2.27 

1.091 

27.68 

0.0550 

0.0221 

50 

2.63 

1.264 

119.5 

0.0638 

0.0950 

75 

3.45 

1.659 

136.6 

0.0835 

0.1083 

100 

3.84 

1.846 

163.6 

0.0930 

0.1300 

150 

5.26 

2.529 

179.2 

0.1275 

0.1472 

200 

6.25 

3.005 

201.1 

0.1515 

0.1600 


























TABLE 4. (Continued) 


Run No. 2 

V m = 41.0 ft/sec (12.5 m/sec) Elevation = 1.848 in. (4.70 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

O 

> 

> 

K f 

Vf/V m 


All 

4.81 

1.000 

0.00 

mn 

0.0000 

10 

4.90 

1.019 

12.80 


0.0102 

100 

5.3 . 

1.102 

118.6 

■nm 

0.0945 

200 

5.8 

1.206 

216.7 

0.141 

0.172 

300 

7.2 

1.497 

261.8 

0.175 

0.208 

400 

8.6 

1.788 

292.2 

0.210 

0.232 

500 

9.8 

2.037 

320.6 

0.238 

0.255 

600 

12.2 

2.536 

309.0 

0.336 

0.246 


Run No. 3 

V m = 35.6 ft/sec (10.8 m/sec) Elevation = 0.898 in. (2.28 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

O 

> 

> 

K f 

Vf/V m 


All 

10.0 

ID 

0.000 

m 

0.0000 

10 

10.0 

BESS 

6.280 

Wmm 

0.0050 

100 

10.0 


62.83 

mmm 

0.0500 


10.6 

mmm 

118.6 


0.0943 

300 

11.6 


162.5 

'if 

0.1294 

400 

13.1 

mmm 

191.8 

0.376 

0.1527 

500 

13.9 

1 .390 

227.0 

0.390 

0.1810 

600 

15.5 

1.550 

243.2 

0.435 

0.1935 

700 

16.7 

1.670 

263.3 

0.469 

0.209 

800 

17.5 

1.750 

287.2 

0.491 

0.228 
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TABLE 4. (Continued) 


Run No. 4 r • 

V m = 28.4 ft/sec (8.66 m/sec) Elevation = 0.388 in- (0.985 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

O 

> 

> 

Kf 

■ v f/ v m ; 


All 

• 9.5 

1.00 

' 0.00 

mm 

0.0000 

10 

9.5 

1.00 

6.55 

mm 

0.0052 

100 

10.0 


62.83 

Wmm 

0.0500 

200 

10.4 

1 HI 

120.83 


0.0955 

300 

11.5 

1.211 

163.9 

0.405 

0.1297 

400 

12.6 

1.326 

199.5 

0.444 

0,1590 

500 

13.3 

1.400 

236.2 

0.468 

0.1 8 80 

600 

14.7 

1.547 

249.7 

0.517 

0.1980 

700 

15.5 

1.632 

283.8 

0.546 

0.225 ■ 


Run No. 5 

V m = 24.0 ft/sec (7.33 m/sec) Elevation = 0.1 58 in. (0.402 cm) 


Convective 
Frequency : 

Convective 
Velocity, ' 
fps 

(m/sec X 3.28) 

O 

> 

> 

K f 

V V m ^ 


All 

9.6 

1.00 

0.00 

0.400 

0.000 

10 

9.6 

1.00 

6.55 

0.400 

0.0052 

100 

9.6 

1.00 

65.45 

0.400 

0.0520 


9.7 

1.010 

129.55 

0.404 

0.1030 

fa* El. 

10!4 


181.25 

0.433 

0.1440 

400 

11.5 

1.198 

218.6 

0.479 

0.1735 

500 

12.75 

1.328 

246.4 

0.531 

0.1970 

600 

13.5 

1.406 

279.3 

0.562 

0.222 

700 

14.3 

1.490 

307.6 

0.595 

0.244 

800 

14.6 

1.521 

344.3 

0.608 

0.274 1 

900 

15.9 

1.656 

395.2 

0.662 

0.315 
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TABLE 4. (Concluded) 


Run No. 6 

V m = 20.0 ft/sec (6.1 m/sec) Elevation = 0.033 in. (0.0839 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

O 

> 

> 

K f 

V f/ V m 


All 

8.3 

1.00 

0.00 

mm 

0.0000 

10 

8.3 

1.00 

7.57 


0.0060 

100 

8.3 

1.00 

75.7 

P i 

0.0602 


8.6 

1.04 

146.1 


0.1163 

300 

9.3 . 

mEM 

202.7 


0.161 

400 

10.4 

UMfl 

241.7 

BE 

0:192 

500 

11.6 

1.398 

370.8 

■sill 

0.215 

600 

12.8 

1.542 

294.5 

0.640 

0.234 

700 

13.1 

1.578 

335.7 

0.655 

0.267 

800 

13.5 

1.627 

372.3 

0.675 

0.296 

1000 

15.0 

1.807 

418.7 

0.750 

0.332 


Run No. 7 

V m = 15.4 ft/sec (4.7 m/sec) Elevation = 0.008 in. (0.0215 cm) 


Convective 

Frequency 

Convective 

Velocity, 

fps 

(m/sec X 3.28) 

o 

> 

> 

Kf 

£ 

> 

> 


All 

5.9 


0.00 

mm 


10 

5.9 

1 

10.65 



100 

5.9 


106.5 

0.383 


200 

6.3 


199.5 

0.409 

1 . 1 I H 

300 

6.3 

1.068 

299.2 

0.409 

0.238 

400 

6.7 

1.136 

375.1 

0.435 

0.298 

500 

7.1 

1.203 

442.5 

0.461 

0.352 
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Figure 24. Convective velocities. 


TABLE 5. CONVECTIVE VELOCITIES; X= 18 ft (5.5 m) 
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Figure 25. Correlation measurement of convective velocities as a function of frequency. 



OUTER PORTION OF BOUNDARY LAYER | LAYER 


o 
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Figure 26. Ratio of convective velocity to mean velocity versus number of oscillations between sensors. 


TEST 1 
NEAR WALL 
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Figure 27. Composite of convective velocity versus number of oscillations between sensors. 


Figure 27 is a composite of the four tests, the results of which are far from conclusive. 
Another possible explanation for the convective velocity being a function of spatial 
separation is that, for closely spaced sensors, a random walk phenomenon is occurring, and a 
conglomeration of many velocities is being measured. With longer and longer separations, 
only the trajectories moving directly in the direction of the two sensors would be correlated. 
The use of different separation distances might be a tool for measuring selected particle 
speed. Other authors have identified this problem [49, 57] . Before a firm explanation can 
be given, much more research needs to be done on this problem. 

The convective velocity below the production zone was examined as a function of 
downstream distance. Figure 28 gives tabular and graphical results for the region near the 
surface. The model is consistent along the X-direction both above and below the production 
zone; that is, the smaller scale structures more nearly approach the mean velocity than do 
larger scale structures. The filtering of convective velocities is as predicted by the model. 


D. Limiting Convective Velocities 

The purpose of this section is to provide experimental evidence to support the 
postulate in Section III that the convective velocity must approach the mean velocity 
throughout the boundary layer after a sufficiently long dwell time. 

Example l of Section III predicted that the convective velocity will approach the 
mean velocity for a boundary layer that had developed for approximately 28 ft (8.55 m) 
with a free-stream velocity of 20 ft/sec (6.1 m/sec). Since the meteorological wind tunnel is 
capable of producing the large growth lengths necessary to test the prediction, the 
convective velocity profiles were measured for a boundary layer developing for 30 ft (9.15 
m) in length. A free-stream velocity of about 25 ft/sec (7.63 m/sec) was used for the study. 
Figure 29 shows the convective and mean velocity profiles at the 30-ft (9.15-m) station. The 
variation of convective velocity, as a function of frequency, was very small. A typical set of 
correlation curves is shown in Figure 30. Figure 30 indicates that all the frequencies are 
moving at very nearly the same velocity. According to Appendix D, paragraph F, one would 
expect a nearly symmetrical probability density distribution for the velocity fluctuations. 
This in turn would produce a symmetrical cross-correlation. Figure 30 shows very 
symmetrical correlation curves at all frequencies. 

Since the formation zone exists along the entire boundary floor, there will be some 
packets being continually fed into the flow all along the boundary layer. This means that 
there will be a few packets that are somewhat newer to the flow than others. Thus, 
complete conformity between the convective velocity and the mean velocity may not ever 
be exactly obtained. 

The present measurements again confirm the physical model. Table 6 gives the 
measurements of Figure 29 in tabular form. 

Very near the floor, where the production zone is still active, one would still expect 
to see some difference between the convective arid the mean velocity even at the large 
X-distances. The region near the floor at a test location 18 ft (5.5 m) downstream shows 
that the convective velocities near the floor are still very much different from the local mean 
velocity (Fig. 28). . 
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ALL MEASUREMENTS TAKEN VERY NEAR THE WALL 

BRACKETED NO. ON CONVECTIVE VELOCITIES DENOTES HIGH PASS FILTER SETTING IN Hi 

NUMBERS ABOVE DATA POINTS INDICATES CONVECTIVE. VELOCITIES OF DATA WITH 
HIGH PASS FILTER SETTINGS OF THAT NUMBER 
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Figure 28. Convective velocities near the surface. 



Figure 29. Convective velocities 30 ft (9.15 m) from roughness. 

The limiting case of the convective velocity catching up with the mean velocity may' 
give some insight into why similarity solutions do not exist except for very thick boundary 
layers. 


E. Convective Velocity Measurements from the Phase Spectrum 

This section gives a brief account of experimental data taken from the phase portion 
of the spectrum to determine transit times and convective velocities. 

Appendix D, paragraph C, gives the theory which allows the phase spectrum to be 
used to evaluate convective velocities and transit times. Equation (D-50) shows the relation 
of transit time between sensors to the shift in phase. Figure 3 1 shows the phase spectrum and 
the calculation of transit time from the phase spectrum, 6 = 2nf t . The data are also given 
in tabular form in Figure 3 1 . Equation (D-5 1 ) gives the necessary relation for determining . 
the convective velocity from the phase spectrum, 
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TABLE 6. CONVECTIVE VELOCITIES; X = 30 ft (9.1 5 m) 
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Figur.e 32 gives a tabular and graphical result of one of the test runs. 

This method is valuable for measuring the convective velocity because one can 
determine the bandwidth of the spectrum for which convective velocity measurements are 
meaningful. That is, only when the time delay determined by the phase spectrum is a 
smooth continuous function would the convective velocity measurement be significant. The 
graphical display of the phase spectrum, phase convective velocity, and phase time delay 
appears, to be. a useful and easy tool for .determining, the significance of a measurement for 
any particular frequency bandwidth. 



0 200 400 600 800 1000 

FREQUENCY, f ( Hz) 



Figure 32. Phase velocity. 
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F. An Exploratory Test 


Convective velocity profiles were taken at a station 1 ft (0.305 m) downstream of 
the roughness near the meteorological wind tunnel inlet. It was felt that inlet accelerations 
would have an effect at this location. The mean, free-stream velocity in the test section was 
approximately 38 ft/sec (11.6 m/sec). Figure 33 shows the results of the measurements in 
the boundary layer. Figure 34 shows the correlation curves over the frequency range of 
interest for one location., 



Figure 33. Convective velocities at entrance of meteorological wind tunnel. 

. j ■ t 

It is seen that the convective velocities are larger than the local mean velocity. The 
higher the frequency is, the more nearly the convective velocity approaches the local mean 
velocity, as predicted by the model. 

This test is presented as an exploratory case where not all of the phenomenological 
aspects of the inlet accelerations are known. 
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Figure 34. Correlation measurement of convective velocities at entrance of meteorological wind tunnel. 



G. Some Thoughts Concerning Eddy Viscosity 

Taylor stated, “It seems natural to suppose that eddies will transfer not only the 
heat and water vapour, but also the momentum of the layer in which they originated to the 
layer with which they mix” [17] . Similarly, it seems logical to describe a viscosity created 
by turbulent interactions by a velocity profile of the turbulent motion rather than the mean 
motion. One might think of this as macroscale interactions of turbulent packets, similar to 
the concept of dynamic viscosity of a gas being caused by collisions of molecules. Describing 
an eddy viscosity by 



where e is the eddy viscosity, seems realistic. Rotta notes that for nonequilibrium 
boundary layers, the concept of a constant eddy viscosity defined by 




becomes more realistic with distance downstream [72]. The results would be consistent 
with a constant eddy viscosity based on the convective velocity profile, because farther 
downstream the convective velocity profile more nearly approaches the mean velocity 
profile. At the distance downstream where the two profiles assume the same shapes, a 
constant eddy viscosity concept for the outer portions of the layer would be identical for 
both convective and mean velocity profiles. Rotta also points out that, for more accurate 
results, the effect of upstream history on the turbulent motion must be taken into account 
[72]. Both the mathematical model and the experimental results show that the convective 
velocity depends upon the upstream history. Since the mean velocity profile may not be 
representative of the turbulent interactions, the concept of eddy viscosity being defined by 
the convective velocity profile of the turbulence is presented as an alternative approach. 


VII. CO NCL US 10 ISIS 


A physical model is proposed to explain the formation and motion of packets of 
vorticity within the turbulent boundary layer. It is hypothesized that a production region 
exists where turbulent packets are . formed. These packets migrate into that portion of the 
flow field which is above and into that portion of the flow which is below this formation 
zone. A mathematical ,/nodel describing the above action was used to predict the convective 
velocity behavior across the turbulent boundary layer. Experiments were then devised to 
check certain of the conclusions derived from the model. 
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The three hypotheses of the model may be summarized as follows: 

1. The undulating outer edge of the boundary layer is the driving function for the 
flow regime near the wall. 

2. Packets of vorticity are produced in the. region which correlate with the outer 
edge variation. 

3. The packets migrate into the boundary layer with some initial vertical velocity, 
their horizontal velocity being that of the zone from which they are produced, 

For the third hypothesis the mathematical model predicts the following: 

1. Convective velocities outside the formation zone are different from the pointwise 
mean velocities. 

2. Convective velocities are filtered by the boundary layer such that the higher 
frequencies are initially accelerated at a greater rate than the lower frequencies. Thus, 
initially, the boundary layer separates the various frequencies, allowing the higher 
frequencies to more closely approach the pointwise mean motion in the boundary layer. 
Convective velocities are functions of frequency. 

3. The acceleration of the packet is a function of the square of the difference 
between its velocity and the velocity of the mean motion of the surrounding fluid. 

4. The limiting boundary layer will have a growth that is linear with respect to the 
stream wise direction. 

5. For fluid structures of all sizes, the convective velocity will approach the mean 
velocity profile for very long growth distances of boundary layers. This is because the mean 
velocity is the driving function on the fluid structures and, thus, will accelerate the packets 
and will be the asymptotic value for the convective velocities. 

The conclusions of the experimental portion of the study may be summarized as 
follows: 


1. There is a correlation between the velocity fluctuations of the highly intermittent 
outer edge of the boundary layer and those very close to the wall. 

2. A production zone exists where fluid packets appear to be formed. Here the 
convective velocity is independent of wave number. The convective velocity equals the 
pointwise mean velocity for this region. 

3. The convective velocity becomes a function of wave number outside the 
production zone. That is, the boundary layer appears to filter the various sizes. The smaller 
sizes are accelerated more than the large sizes. This may also be stated as follows: packets 
whose characteristic frequency are higher are convected at a rate more nearly the mean 
velocity than are packets with lower characteristic frequencies. 
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4. For boundary layers which have developed over long lengths, the convective 
velocity becomes independent of wave number in the outer regions due to the asymptotic 
relationship between the convective and mean velocities in the outer regions of the 
boundary layer. 

5. The probability density distributions of the turbulent fluctuations are skewed 
when the convective velocity does not equal the local mean velocity. 

6. The cross-phase spectrum may be used to determine the uniformity of 
convection velocity over wide bandwidths. 

The experimental results are consistent with the requirements of the physical model. 
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APPENDIX A 


SOLUTION OF EQUATIONS (10) AND (11) 

The equation needed to solve the X-direction of motion of the packet of vorticity is 
equation (10). That is, 

dv x c d y 2c d v m v x c d v m > 

dt X i X i X i 

Noting that Cj , , and V m may all be functions of time, equation (10) becomes the 

generalized Riccati equation. This is a special case of the Abel equation [63, 64], The 
Riccati equation is a differential equation of the form 


y(t)' = h(t) + g(t)y + f(t)y n . (A-l) 

Kamke [63] states that the solution may be found, if suitably selected 
constants a and P may be found, where 


(w)^ = IP + af h(t)e' / g (t)dt dt] 


(A-2) 


such that Z(t) = 



is a solution of the linear differential equation, 


Z(ty - g(t) Z(t) = ah(t) 


(A-3) 


Then, the solution of the original differential equation follows from 



(A-4) 
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where U(t) is determined by 



(A-5) 


For equation (10), the constants a and P must be selected such that 


m 


= e 


Xi 


[0 + a / 


r ,V 2 
^d v m 

X, 


;^d, 

a A. i 


dt] 


(A-6) 


Equation (A-6) requires an explicit relationship for , V m , and X i as functions 

of time. For many cases, these will not be known and approximations must be used. One 
may wish to consider , and V m to be independent of time. The differential 

equation then becomes a special form of Abel’s equation to which the solution is given 
below. 


By equation (A-l), 


. . c dV 

h - -TT - 


; g = 


- 2C d v m 

X, 


; f = 


x, 


and 


- w - 


m 


Using (A-2) 


V = e J X! 
v m e 


/SV dt 


[p + 


r c dV /-4r^ dt , tl 

a J — r e dt] 

A i 


(A-7) 


(A-8) 


(A-9) 
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At t = 0 , V m = V m => 0 = V m . Differentiating equation (A-9) yields 


0 = 


- 2 C d v m 2 

X, 


“ c dV 

h 



' . (A-10) 


From (A-5), 


U' 


% 

Xi 


(U 2 - 2U + 1) 



and 


dU 

U 2 -2U+1 


1 

1 -U 


V m C d 


t + C 
1 


(A-ll) 


Equation (A-4) gives V x = ZU , where V x = Vj at t = 0. The equation for velocity of 
the packet besomes 


V x = 


m 


x »( V m- V i) 
tC d( y m- V i) + 


(A-12) 


Integration of equation (A-12) yields the X-direction of motion for the packet: 


X,(V m -Vj) 

X(t) = V m t - ” _ V .) ' gn[Xl + tG d( y m - y i>] + C • .(A-13) 

The initial condition that X = 0 , at time t = 0 , yields 

X(t) = V m t - -^8n[X, + tC d (V m - Vj) + . . . (A- 14) 
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Thus, the equations of motion for the packet in. the X^direction are as follows: 
Acceleration 


d 2 X(t) 
dt 2 


dV x (t) 

dt 



2 V V 
^ v m v x 


+ V 2 1 
m 1 


(A-15) 


Velocity 


M = v {t ) 

dt v xt l) 


m 


^( v m- v i) 

x. + tc d (v m -v i ) 


(A- 16) 


Distance 


X(t) = V m t - 8n[X, + tC d (V m - Vj)] + . (A- 17) 


The technique used to solve this equation was selected to give reference for general 
use. The following is given as proof that the previous approach will give a solution for the 
equation : 


Y' = f(x) Y n + g(x) Y + h(x) 


(A- 18) 


where 



e [p + a f he 'f&x dx] 


(A- 19) 


and 


r i « + c 

J U n - «U + 1 


* /(■ f) /n hdx 


(A-20) 
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The hypothesis being that Y = ZU is a solution of the differential equation, (A-l 8). 

Proof 

Z' = e^S^ x g Ze'-^ x + e-^ x ahe“-^ x = gZ + ah (A-21) 

U ' ■ T< U " - “ U + *> ■ (A-22) 

From the hypothesis Y = ZU, 

Y' = Z'U + U'Z . • . ■ ■. ’! : (A-23) 

Equation (A-23), combined with equations (A-21) and (A-22), yields: 

Y' = U(gZ + ah) + Z (U n - aU + 1) 

Y' = gZU + hU n + h ' ■ ' 

Y' = hY n + gY + h (A- 24) 

• , * ; * . • . •> ■ ■ 

Equation (A-24) is the original equation, and thus Y = ZU is a solution to equation 
(A-l 8).* 

The solution to equations (10) and (11), where , V m , and Ai are independent 
of time, may also be found by noting that equation (10) may be written as 

dV„ 

’ - C(V X 2 + 2aV x + b) = 0 , (A-25) 

*Proof of existence thanks to Dr. H. G. Krause, NASA, MSFC. 
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where 



-V m ; and b 


V 2 
v m 


Equation (A-25) is now a special form of Abel’s equation [63] where the solution is 

known 




+ C) 


(A-26) 


Applying the initial conditions V x = Vj at time, t = 0 , 


V - V = 
m v x 


^ V m- V i> 
tC d( V m-Vi) + 


(A-27) 


which is the same as equation (A-l 2) derived earlier. 

Equation (A-27) may now be viewed as a differential equation to describe the 
distance the packet moves in the X-direction. 


dx *.( V m- V i> 

dt -tC d ( V m-Vi) 


(A-28) 


the solution being 


x 


-^( V m~ V i) 

C d ( V m- V i) 


Sn[X t + tC d (V m - Vj)] + V m t + C 


(A-29) 


which is the solution to equation (11). 
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Applying the initial condition that X = 0 ‘at time t = 0 , we obtain 


x(t) 


= -£-en[x, + tC d (V m - Vj)i + v m t + 


(A-30) 


equation (A-30) being identical to equation (A-17). 

For easy reference the basic equations for the X-direction of motion of the packet 
are given. (Note: , Xi , and V m are independent of time.) 

Acceleration of the packet in the X-direction 


A x 


dVx Cd(Vm-Vj) 2 

dt X, 



- 2V V 
i4V m v x 


V 2 1 
v m 1 


(A-3 1 ) 


Velocity of the packet in the X-direction 


V = V 
v x m 


X ^ V m- V i) 

X, +tc d (v m -v i ) 


(A-3 2) 


Displacement of the packet in the X-direction 


x(t) = V m t --^i-gntX, + tC d (V m - V t )] + - X ^ .L (A-3 3) 
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APPENDIXB 


SOLUTION OF EQUATION (28) 

l 

The solution of equation (28), 


_ W 

dt 


2C, 


^xf-ICvV + VilVx +-XT-KVV + Vi ] 2 ■ 



- (CyVyjt + Vj )] 2 



is performed by using a change of variables. Let 


Z = [V x - (C v v yi t + Vj)] 


(B-l) 


Then, 


dZ 

dt 


dV x 

dt 


C V 
v 'v v yi 


(B-2) 


Equation (28) may now be written as 


dZ = c d z2 r r 

dt \i ~ v yi 


(B-3) 


Solutions to this differential equation may be found in Reference 63, p. 296, no. 
123, or Reference 64, p. 229, no. 55. 

Murphy [64] gives for y' = a + by 2 , 

(i) ry = a tan (C + rx) ; r = \/aF 
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(ii) sy = a tanh (C + sx) ; s = v^ab 

(iii) y(C - bx) = 1 ; a = 0 

(iv) y = ax + b;b = 0 

For the case in question, y = z, x = t , a = -C v V y j , b = C d /Xi , 
and ab = (-C d C v V yi )/Xi . > 

1 

The velocity profile in question has velocity increasing with elevation; thus, C y > 0 . 
Also, Vyj > 0 , C d > 0,and\i > 0 . Equation (ii) must therefore be used. 

sy = a tanh (C + sx) ; s = V-ab 


Therefore, 


^ c dC»V y i 




( V x * 


C V t 
W v yi l 


- vp 


-C v V yi tanh 



Solving for V x , we obtain 


-C v V yi tanh 

c + ( 

_ 

^dCyVA 
^ X, ) 

\ h 

' t 

+ C y Vyt + Vj 

( 

V*v y i) 

V ) 

r 


At t = 0 , V x = Vj =»• C = 0 . Thus, the equation for the velocity in the X-direction 
becomes 


V v = 


_c V • 
-v- v v yi 


(c d c v v yi /x,) 


Vi 


tanh 


CjC V • 
v d v yi 

X, 


\‘/2 


CyVyit + 


v i 


or 


(B-4) 
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V = V 
v x m 


C V • 
• v v yi 


(c d c v v yi /x,) 


, - ■ 


Integration of equation (B-4) yields 


x(t) = C y V yi -j- + Vjt --^-8n cosh I 


CaCL'V. 


iV2 


*501 t + c . 


At time t = 0 , x(0) = 0 =» C = 0 , which makes the X-direction of motion become 


, X, „ „ l C d C v V yi f , 

m t - fT-Cn cosh 1 - M t 


x(t) = V 


X, 


(B-5) 


Since V m = C v V yi t + Vj , C v V yi = (V m - Vj)/t . 

The equations for velocity and distance of the packet . may now be written as 
follows: 

Velocity of the Packet 



Translation of the Packet 


x(t) = 


V m + V i + . X, 


t - -pr- £n cosh 


_SL(v _ y.) t 

TP 1 m 1 


V2 


(B-7) 
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APPENDIX C 


SOLUTION OF EQUATION (10) 

Equation (10) for a velocity gradient of V-^ = (C y y) ' + Vj , where y(t) is 

described by V y jt = y(t) , will be solved by the Kamke approach outlined hi Appendix A. 

For the Kamke approach, 


h = 


C ,V 2 
^d v m 

Xi 


~ 2 Cd v m 

Ai 


; f = if ; and z = (t) 7 ' 


= V 


m 


Given the velocity gradient, 


m 


= (C v V yi t) 


V, 


+ v i 


Z' becomes 


(C-l) 


Z' = -i-(C v V yi t)' 6/7 C v V yi = y (v m ~ v i) 6 c v v yi • ( C ‘ 2 ) 

However Z' also equals gZ + ah ; thus 


Z ' = T (V m- V i)’ 6C v V yi = 


-2CjV 2 

z ^d v m 

Xi 


+.. a 


C ,V 2 
^d v m 

'Xi 


(C-3) 


Solving for a , we obtain 
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« = T< V m - v i'‘‘ c v v vi + 


2C d VJ 


m 


yi \i 


CjV 2 
^d v m 


■ J 7 ["■. v ^i"'] 
* 


% 2Cj r v 

C„v„ ; + [5 C v V vi t} 7 


-V yi X, I v y 1 


x, 


C d[(C v V yi t) + Vj]' 


At t = 0 , 


a — +2 ; a 2 - 4 = 0 
Solving for U , we obtain 


/• 


dU 

U 2 - aU + 1 


+ C 


r dU 
J (1 -U) 2 


+ c 


1 

1 -u 


+ c 


From equation (A-5) 


1 

1 -U 


+ 


C 


- TT /t(c v v yi t) v ’ 


+ Vj] dt 


Integrating and solving for U ; 


U 


C d r 7 

X, L8 



C V -t 

'- V v yi l 


% 


+ Vjt - c 


Using equation (A-4) 



CC-4) 
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Using the initial conditions that V x = Vj at t = 0, equation (C-4) may be solved 
for the constant c yielding 


c 


m 


(Vm-Vi) 


Equation (C-4) for the velocity in the X-direction becomes 


V = V 

x v m c 


v m< v m- v i> 


if < V m - V i> [l C » V y.> % + Vjt] + V m 


(C-5) 
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APPENDIX D , 

SIGNAL ANALYSIS 

* 

A. General Analysis 


The objective of this section is to establish the analysis which will be used, to extract 
information from time series of finite duration. For convenience, it may be supposed that 
the signal is coming from a hot wire anemometer, since this is the type of signal used in the 
experimental portion of this paper. The figure below is given for illustration. 



The recorded signal is of some given finite length, (T 3 , T 2 ). The present examination 
will concern itself with less than the total record length; for example, 0 < T < T^ 

Common statistical parameters that may be used to analyze such a function are 
generally based on the first four moments of the function: 

1. Mean value, 


M f(t) 


f(t) 



f(t)dt 


2. rms value, 


ff f(t) 



(f(t) - M f(t) )’ dt 


'/2 


3. Skewness factor, S^: 


(D-l) 


(D-2) 
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T 

T, { [ f{*> - M fct) l> , 

Sf(t) = J fT 1 (D-3) 

T, { ( f <» - dt J ! 

4. Kurtosis factor, k 


T 

(f(t) - m to )* it 

k f«o * r- — t; T7 • 

T7 { («*> - “fd)) 1 dt J 

5. Auto-correlation, Rff(r): 


%f(r) 


T 

4- /' f(t)f(t + r)dt 
11 0 


(D-5) 


6. Auto-covariance, c^t): 


c d r)= 'T l f l tf(t) - M f(t) ] [f(t + r) - M f(t+T) ]dt . (D-6) 


It is easily shown that only one of the mean values need be used to calculate the 
auto-covariance function; that is, one of the mean values may be set at zero or any 
convenient constant value while using the true mean for the other mean value. 

7. Normalized auto-covariance, Cffir): 


„ , , _ c ff< 7 > 

Cff( T ) — - 

ff f(t) °f(t) 


(D-7) 
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The above expressions give meaningful statistical values in either the time domain or 
time-delay domain. It will also be important to know certain statistical values of the original 
function in the frequency domain. Taylor [ 13] -elucidates on work by Simmons in 1938 to 
show the formal beginning of the importance of the frequency domain. 

The simplest approach for basic understanding is to examine the function f(t) and 
find the best-fit sine wave with n periods in time Ti . A least-square fit from 0 to T, , 
for n oscillations gives 


dt = 0 . (D-8) 

Using the Leibniz theorem for differentiation of an integral in equation (D-8) yields 


_d 

dB 


/' 
n n 


B n sin _ f( t ) 



2B n sin 2 


27rnt 

T, 


2 sin 


27rnt 

Ti 



dt 


0 


(D-9) 


Integration of equation (D-9) gives 


T, 

B„T, - 2 f 
0 


sin ?jr— • f(t) dt 
t i 


= 0 


(D-10) 


Solving for B n , we obtain 


D — * 

n t,/2 


J\ sin f(t) dt 


.0 


Tx 


(D-ll) 


B n is now recognized as being the nth Fourier coefficient of the sine terms of -f(t) . The 

above procedure may also be used for all cosine terms with integer oscillations over the 
interval (0,Tj). 
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A more satisfying approach is to assume that f(t) over (0,T, ) may be represented 
by a series of harmonic components, such as 


f(t) = fj 


i=0 


A n cos 


27rnt 

T, 


+ 


B n sin 


2nnt 

T, 


(D-12) 


To solve for any A n or B n , multiply both sides of the equation by the sine or 
cosine function which accompanies the A^ or B n , and then integrate both sides over the 
interval (0/1^ ). The integration is solved, leaving 



cos 


27rnt 

T, 


f(t) dt 


and 


(D-13) 


B n T, /2 



sin 


27rnt 

T, 


f(t) dt 


(D-14) 


Using the above equations leads the right-hand side of equation (D-12) to converge 
absolutely and uniformly to f(t) over the prescribed interval if f(t) is continuous 
and f'(t) is at least sectionally continuous over the interval (0,T, ) . The only exception is 
at the end points where the right-hand side of equation (D-12) yields y2[f(0) + f(T x )] . 
The above restrictions are applicable since the functions that are dealt with are both 
continuous and continuously differentiable. Thus, the Fourier series will converge within 
any given error with a finite number of terms. 

Equation (D-12) may be written in the form: 


f(t) 



i27rnt 


C n e 


T, 


(D-15) 
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where 


C n = 


1 T > 
11 0 


-i2?rnt 


f(t)e 


T, 


dt 


Noting 


u = 


2nn 

T, 


coTi 

n = 


2n 


and Aco = [(n + 1) - n] = 

ti 1 1 


equation (D-16) becomes 


f(t) 




Acj 


If the time period becomes long, Ao> approaches du> with 


- / 


T|C n 

2jt 


iwt 


dco = f g(co )e 1W 1 do> 


iwt 


and 


T, 


gfw.) = T t C n = f f(t)e' lwt dt 
0 


Cross-channel manipulations require the following statistical concepts: 
1 . Cross correlation of f(t) and g(t) : 


(D-16) 


(D- 17) 


(D-18) 
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(D-19) 


R f g ( t) = ~ / f(t)g(t + r)dt 


2. Covariance of f(t) and g(t) : 


c f g ( T > 


T 

f\ [f(t) - M f(t) ] [g(t + r) - M g (t + r)].dt 


(D-20) 


3. Normalized covariance of f(t) and g(t) : 


Ce jr) - -!fig 

O f(t )O g (t+T) 


(D-21) 


. . It is permissible that the functions of f(t) and g(t) be moments of functions, such 
as f(t) = h(t) n . 


B. Restrictions Caused by Physical Limitations 

The previous section on signal analysis provides the mathematical tools with which 
to analyze real time functions. This section shows the feasibility of defining turbulent flow 
structure with the preceding statistical quantities in view of physical restrictions placed on 
the flow. It is also an objective of this section to show how restrictions placed on a single 
time function propagate the function’s restrictions into moments and cross-moments with 
other functions. 

For the purpose of analysis, consideration will be placed on a finite flow field with 
finite velocities. 

Consider two points in the finite flow field domain: A at (x, and B at 

(x 2 ,y 2 )- Consider also the velocity of the fluid at point A in the direction a and the 

velocity of the fluid at point B in the direction b. For simplicity, it is assumed that the mean 
flow is in the direction of the X-axis with a constant mean value. Then, 


u a (t) = U x e ax + u a (t) at A 


(D-22) 
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and 


U b (t) = U x e bx + u b (t) at B , (D-23) 

where the e^’s are the direction cosines of the ith direction with the jth direction. 

The cross-correlation becomes 

R u a u b (') = U a (t) U b (t + r) = U x 2 e ax e bx + u a (t) u b (t + r) (D-24) 


If the signals are ac-coupled, that is, if only the fluctuating portion of the signal is of 
concern, then, 


R u a u b ~ u a^ u b^ + r ) • 


(D-25) 


Since u a (t) is the fluctuation portion only, i.e., u a (t) = 0, the Fourier 
representation of u a (t) over the interval (0,T i ) may be written as 

u a (t) = E ^ a n c °s -^jr— ^ + b n sin^^-^ , (D-26) 

where 


®n 


2 

T, 



, , , 2jrnt j. 
u (t) cos -7= — -dt 
m 


(D-27) 


and 


105 



b„ = 


■ 4-7' 


2nnt 


J u b (t)sin dt 


(D-28) 


Similarly, 


%(t + t) 


- 1 k 

n=l \ 


cos 


27rn(t + r) 

T, 


b' n sin 


27rn(t + r) \ 

T. ) 


(D-29) 


with 


2 Tl+T 

n “ \ I u b(t + r) 


cos 


2an(t + r) 

t; 


dt 


(D-30) 


and 


b' = 


n 


T, +T 

2 r \ • 2im(t + t) 

Y J u b (t + t) sm — -dt 


(D-31) 


Letting t' = t - r for any r , dt' equals dt , 


T, 


t Z r /,K 2jrnt ,, 

n ~ ~T J Ub ^) 005 dt > 


T, 


T, 


(D-32) 


and 


T 

u> 2 r' • 2nnt' A . 

b n = X J ~tT dt 


(D-33) 
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For any given r 


^(t) u b (t + r) 





, 2unX . 
a _ sin — ■ sin 

*i 


27rnt 

T, 


+ 


b' 


sin 


2?mt 

T, 


cos 


27rnt 

T, 


. , 27mt • 
b n cos sin 



dt 


(D-34) 


Orthogonality of the function is noted and the integration is performed: 


u a (t) u b (t + r) 


_1_ 

Ti 



a n a n 


cos 


27inr 

Ti 


+ 


II 

2 


a n b' n sin 


27rnr 

T, 


!> 

2 


b n a 'n 


sin 


27rnr 

T, 


+ 


T, 

2 


b n b 'n cos 



1 

2 



a^^ cos 


271 nr 
T, 


+ 


a n b 'n sin 


2-nnr 

Ti 


' b n a 'n sin 


277 nr 
T, 


b n b' n cos 


277 nr \ 

T. } 


(D-35) 
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The average value of the correlation over the r region of interest becomes 


i r m I “ / 

U b^ t + T ^ = 17” I T Z ( a n a n 
m r n-1 \ 


2nm , . u 2nm 

cos — ^ — + a n b n sin — — 


, , 2?rnr , , , , 27rnr \ , 

b n a n sin — ^ — + b b n cos — — 1 dr 


(D-36) 


Integration over symmetrical limits eliminates all odd functions, and integration of even 
functions over symmetrical limits yields twice the value integrated over half limits. 

The average value may now be written (where the presubscript, r , denotes 
averaging over r ) as 


j J m “ / 

T u a (t) u b (t + r) - f Yj ( a n a n 

m o n=l \ 


, . : , , 2rrnr \ . 

+ b n b n cos -j — 1 dr 


(D-37) 


Letting (a n a^ + b n b^ ) = 2 a^, equation (D-37) becomes 


r u a(t) u b (t + r) = y-'f™ Y \ a " n C0S dt 

m o n=l 1 


(D-38): 


~ a" n t 2irnr m 

r u a(t) u b (t + r) = £ Inn sin ~T “ 

n=l m 


(D-39) 
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Each term of equation (D-39) is identically zero whenever 



m = 1, 2, 3, 4, ... 


This zero value indicates that for structures defined within a finite time, their correlations 
are also defined within an identical finite time. 

In equation (D-35), by letting 


(a n a' n + b n b' n ) = 2a" n and (a n b' n - b n a' n ) = 2b" 
the cross-covariance becomes 


^u a u^ T ) U b^t + T ) ^a n cos x, 

+ b" n sin-^^ . (D-40) 

Equation (D40) is identical to equation (D-26) except for the notation r and t . Thus, any 
restrictions of a physical nature imposed upon the original functions are propagated 
invariantly to the correlation function. 

If physical limitations allow the time series to be considered stationary 
through i moments, it is now proposed that products of moments will retain the original 
degree of stationarity. 

Consider U a (t) and U^(t + r) stationary functions through their ith and jth 
moments, respectively; that is, the moments of U a (t) and U^(t + r) are independent of 

the ’placement of the interval to be analyzed as long as the length of the interval stays 
constant. Expressed mathematically, using the Leibnitz theorem, 


a+T° 

U a °(t)dt = 0 = U a n (a + T 0 ) - U a n (a) (D-41) 
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and 


- <*+T 0 

f U b m (t + r)dt = 0 = U b m ( a + T 0 +r) - U b m (<* + r) , (D-42) 
a 


where n and m represent any stationary moments of u a (t) and U b (t + r) , respectively. 

Examining the product of stationary time series for stationarity, the Leibnitz 
theorem becomes 


a+ T° 

/ U a n (t) U b m (t + r)dt = U a n (a + T 0 ) U b m (a + T 0 + r) 
a 

I 

- U a n (a) U b m (a + r) 


= U a n (a) U b m (a + r) 


- u a n (a) u b m (a+r) . (D-43) 

A corollary of equation (D-43) could be expressed as follows: “The product of any 
number of stationary functions is also stationary.” Thus, equation (D-43) would be 


a+T 0 

^ / U a n (t)U b m (t) ... UJdt = U a n (a + To)U b m (a + T 0 ) ... UJ(a + ^) 


-U a n (a)U b m (a) ... UJ(a) . (D-44) 
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However from equation (D-41), U^Ca + T 0 ) equals UJ(a) . Equation (D-44) now becomes 

identically equal to zero, and thus the product of moments of stationary functions is also 
stationary. This corollary may be expressed as follows: “The product of moments of two or 
more stationary functions is independent of the placement of the integral interval.” Thus, 
by definition, the products of moments of two or more stationary time functions are 
themselves stationary. 

The above analysis indicates that the product of stationary functions is only 
stationary to the degree that each individual function is stationary. 

j The conclusions of this section may be expressed as follows: 

1 . Any structure or function described in a finite time may also be described by 
correlation functions with finite time delays. Expressly, the time-delay range must be of the 
extent of the structure or function for a complete correlation description. 

2. The average value of the correlation function is a function of the correlation’s 
upper time delay and is only zero at prescribed time-delay values, equation (D-39), 
excluding, of course, the trivial case where all Fourier coefficients are zero. 

3. The product of moments of two or more stationary time functions is also 
stationary to the degree of the least stationary component. That is to say, the degree of 
stationarity of a product of stationary functions is the same as that of the least stationary of 
its components. If any moment of any component is not stationary, this product will also 
be nonstationary. 

4. A structure defined in a length of time T, can only be described fully in the 
time-delay domain if r max - T m j n > T } . 


C. A Physical Interpretation of the Fourier Phase Spectrum 

Harmonic analysis is one tool which may be used to present data in such a way that 
their full meaning may be easier or better understood. Phase shifts of individual harmonics 
may be used to ascertain convective velocities. The phase spectrum is treated 
mathematically in this section, along with its physical interpretation. 

Consider two time functions and their respective Fourier transforms: 


x(t) 


oo 



a n cos 


2nnt 

T 


+ 


b n sin 


27mt 

T 


OO 



(a 2 , 


bV 


i/ 2 , 2nnt 

n cos( 


- e 


<y 


(D-45) 
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and 


y(t) 


oo 



C n cos 


2nnt 

T 


'n 


sin ■ 


2rrnt 

T 


op 





(D-46) 

t 


Where 0 xn and 0 yn represent the phase shift of each harmonic at the 
origin, 0 xn equals tan' 1 b n /c n and 0 y>n equals tan' 1 d n /c n . The correlation of x(t) 
and y(t) becomes 


T r 


i*xy( T ) ... t f 
o 


x(t) y(t + r)dt = 


ap c 0 
4 



a n c n cos 


2?rnT 

T 


a n d n sin 


27rnr , _ 27rnr . , , . 27rnr 

"f b n c n SU1 + b n d n cos T" 


(D-47) 


The correlation may also be thought of as a function, without regard to the 
functions involved in its formation. The Fourier transform of the correlation function could 
be written as 


R X y( r ) 



+ 


OO 



e n cos 


27rnr 

T 


+ 


fnSin 


27rnT 

T 





\ 



(D-48) 


where 0 cn is the phase shift from the origin of the nth harmonic, 0 cn '= tan' 1 f n /e n ■ 
Substituting from equation (D-47) for f n and e n we obtain 
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(D-49) 


,-i Jn 


a„d„ — b„c„ 
-! n n n n 


d r n = tan' 1 = tan' . , , 

c ’ n e n a n c n + b n d n 


= tan' 


£n_^n 
c n a n 

a n c n 


Thus, 


= tan 


tan0 „ n - tan0 v _ 

-i y ? n x.,11 

1 + tanf? y,n tan0 x,n 


= tan' 1 [tan(0 yn -0 xn )] 


e c,n 6 y,n 9 x,n 


(D-50) 


The physical interpretation of equation (D-50) may be stated as follows: “The 
phase of the nth harmonic of the correlation is equal to the difference of the starting phases 
of the nth harmonics of the original time series which were used to calculate the 
correlation.” A special case of this result occurs when x(t) = y(t); thus, 0 y n equals 0 X n . 

This means that 0 cn must be equal to zero from equation (D-50); that is, the phase 

spectrum for the auto-correlation must equal zero for all frequencies. That the phase 
spectrum is equal to zero for the auto-correlation is well known and may also be arrived at 
by considering the evenness of the autocorrelation. > 

The correlation receives its major contribution from the nth harmonic when 


/ 2rrnr a \ . 

" S (-T- - 9 c,n) * 1 


that is, when 


_ 27mr 
c,n j 


2?rfT 


(D-51) 


If all waves travel such that they peak at approximately the same time delay, the 
correlation curve will have a pronounced maximum. The calculation of the convective 
velocity of any harmonic is then the separation distance between the sensors divided 
by e cA !2iti : 



(D-52) 


V„ = 


2«f{- 


c,n 


Equations (D-51) and (D-52) show that the time delay may also be calculated from the 
phase spectrum; i.e., r = 0 c n /27rf . Figure 31 gives the phase spectrum and time-delay 

spectrum as functions of frequency for a given test. 

The conclusions of this section are as follows: 

1 . The phase of a particular harmonic describing the correlation function equals the 
difference between the starting phases of the two functions used to calculate the 
correlation; i.e., 


^c,n 0 y,n ' fl x,n 

2. The optimum time delay for a particular harmonic to convect from one sensor to 
another is given by the equation 


r p = 9 c,n l 2iri • 

3. The convective velocity of particular portions or individual harmonics may be 
calculated from the equation 


V c = 2yrf£/0 

These relations describe the convective phenomena in terms of phase and frequency. 


D. Special Properties of Correlations of Time Functions 


The time functions used to describe physical models must be of finite duration as a 
matter of practicality. Thus, for the following derivations the Fourier series is used rather 
than the more familiar Fourier integral. The time length analyzed in all cases will be of total 
length Tj . The following cases will be examined: 


1 . Auto-correlation of function f(t) . 

2. Cross-correlation of function f(t) with g(t) . 

3. Cross-correlation of a function, f(t) , with its first time derivative. 
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4. Cross-correlation of a function, f(t) , with a second function’s 
derivative. 

5. The derivative of the auto-correlation with respect to t , time delay. 

6. The derivative of the cross-correlation with respect to r , time delay. 

7. The derivative of an auto-correlation with respect to time. 

8. The derivative of a cross-correlation with respect to time. 

9. The second derivative of a cross-correlation with respect to time delay, 

1 . Auto-Correlation of Function f(t) . The auto-correlation will be 
considering the spectral content of the function to be auto-correlated. 
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(D-53) 
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The auto-correlation may now be written as 


, T, T, 

f(t)f(t + r) = ~ / f(t)f(t + r)dt = y: f 


T, 


0 


*♦■£,(*. “tt 

n=l ' 


. • 2irn r 

* B n sm T: 


r|4 + £, ( A 'n- 


2rrn(t + r) 
Ti 


+ B' n sin 


2?rn(t 




dt 


= -L f 

T, J 


0 


Ap 

2 


£ <*■ 


27rnt . n 27rnt 
cos — ^ — + B n sin — y 


r) 


A' 

2 


+ 2 ( A ” 

n=l 


2irnt 2nn r 
cos ■ ^ cos — 


A , 27rnt 2rrnr 
A n sin —y — sin y ~ 


. . 2jrnt _ 2irm , 2irnt . 2jrnr\ 

+ B n sin —y — cos — y — + B n cos —j — sin --- -- j 


dt 


(D-59) 


Upon integration of equation (D-59), the auto-correlation becomes 
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Equation (D-60) may be rewritten as 
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If the function f(t) is stationary , 


A n = A 'n and B n = B 'n 


Thus, for a stationary function, 


< A n B 'n - A 'n B n) 


= 0 


The auto-correlation of f(t) is 
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And, if the function f(t). is stationary, the auto-correlation becomes 
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f(t)f(t + T) = + l 

n=l 


°° A n +B n 2jrnr 

- cos 


T, 


(D-6 1 ) 


(D-62) 


(D-63) 


1 17 ' 



The spectra of f(t) and f(t) f(t + r) are seen to be related because A 0 2 /4 equals 
(A 0 /2) 2 , and (A n 2 +B n 2 )/2 equals one-half of the power spectrum of f(t) at the 

respective frequencies. 

2. Cross-Correlation of Functions f(t) with g(t) . The functions of f(t) , A n , 
and B n are described by equations (D-53), (D-54), and (D-55), respectively. g(t + r) is 
described by equations (D-56), (D-57), and (D-58), where f(t + r) is replaced by g(t + r). 

The cross-correlation may thus be written as 


f(t) g(t + t) 



f(t) g(t + t) dt 


(D-64) 


The correlation in terms of its Fourier components becomes 
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(D-65) 


If the spectrum of f(t) is considered as A n + i B n and the spectrum of g(t + r) is 
considered as A^ + i B^ , where i = y/T means a phase shift of 90 degrees, the spectrum 

of equation (D-65) may be obtained by multiplying the spectrum of g(t + r) by the 
complex conjugate of the spectrum of f(t) . 

3. Cross-Correlation of a Function, f(t) , with its First Time Derivative . 
Differentiating equation (D-53), 
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where A^ and B^, are defined by equations (D-57) and (D-58), respectively. Thus, 
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Upon integration term by term, 
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- n(A n A' n + B n B' n ) sin . (D-68) 

ll / (Concluded) 

If the function f(t) is stationary, A n = A^ and B n = B^ . Equation (D-68) then 
becomes 
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VM at 
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Equation (D-69) evaluated at r = 0 yields 


f(t) = 0 , for r = 0 . (D-70) 

* ot - . . * • ' . 


4. Cross-Correlation of a Function, f(t) , with a Second Function’s First Time 
Derivative . f(t) is defined by equation (D-53), and g(t + r) is defined by equation (D-56). 
deft + 7) 

The term — is defined by equation (D-66), where f(t + r) is replaced by g(t + t). 

Equation (D-68) now may be written to define f(t) + 7 ^ as 

0 t 


ag(t + r) 

vw at 




n 


- B n A' n ) cos 


27inr 

T, 


~ "(AnA'n 


+ B n B' n ) sin 


27inT \ 

T* 7 ’ 


(D-71) 


5 . The Derivative of the Auto-Correlation ■ with Respect to t , Time 
Delay. Differentiating the auto-correlation with respect to time delay yields 
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Note that 
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3f(t + r) 
3 r 


(D-73) 


From equations (D-72) and (D-73), 


= fm af(t + T) 

a r w at 


(D-74) 


Therefore, from equation (D-70), the auto-correlation must have a maximum at r = 0 ; 
i.e., the slope at r = 0 must equal 0. The spectrum of equation (D-74) would follow from 
equation (D-68). 

6. The Derivative of a Cross-Correlation with Respect to r , Time Delay. In 
equation (D-72), replacing f(t + r) by g(t + r) gives 


9R fg T = fm ag(t + r) 

3t K ’ at 


(D-75) 


The spectrum of the above follows from equation (D-71). 

7. The Derivative of an Auto-Correlation with Respect to Time . Taking the time 
derivative of the auto-correlation function gives 


9 Rff(r ) 
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(D-76) 


Since Rff(r) is not a 
however, give 


function 


of time, its time derivative is zero. This does, 
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Therefore the spectrum of f(t + r) would follow from the negative of equation 

(D-68). 9t 

8. The Derivative of a Cross-Correlation with Respect to Time. Replacing f(t + r) 
by g(t + r) in equations (D-76) and (D-77) gives 


3R fg (r) 

at 


and 


f(t) 3g(t + t) = _ 3f(t)_ ( (D-79) 

at at ’ 


The spectrum of g(t + r) follows from the negative of equation (D-71 ). 

o t 

9. The Second Derivative of a Cross-Correlation with Respect to Time Delay, r . 
Using equations (D-72), (D-73), (D-75), and (D-79), the second derivative of the 
cross-correlation with respect to time delay becomes 


3 2R fg(r) 

3r 2 


a 2 f(t) 

3r 2 
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(D-80) 
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To consider the auto-correlation, merely replace g(t + r) by f(t + r) . 

The spectral content of equations (D-80), (D-81), and (D-82) may be expressed as 
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For an auto-correlation of a stationary function, equation (D-83) becomes 
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Notice that equation (D-83) gives the spectral content for equations (D-80), (D-81), and 
(D-82). Notice also that the spectra of equations (D-83) and (D-84) are identical with the 
spectra of equations (D-65) and (D-63) except that each term is multiplied by the square of 
its angular frequency. Equation (D-80) gives the relationship of the correlation of two 
functions to the correlation of the derivatives of the two functions. 


E. The Effect of Space and Frequency Weighting Functions on the 
Normalized Cross-Covariance 

Let f(t) be the fluctuating portion of the signal which first senses the flow 
phenomena. Let g(t) be the fluctuating portion of the signal which then senses the same 
flow phenomena at some later time, the separation between the sensors being £ . 

Assume for the moment that f(t) is changed only due to spacing and frequency 
before it becomes g(t) . Then, 
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where F(£,f) is the above weighting function in space and frequency and r n is the time it 

takes for the nth harmonic to reach the second sensor. The correlation of f(t) 
and g(t) becomes 


T oo 

f(t) g(t + r) = y / 2 (a n cos nut + b n sin nut) 

0 n=l 


oo 

F(S,f) Jj [" a n cos nu (t - t + r) + b n sin nu(t - r n + t) 1 dt 
n=l L J 


°° a 2 + h 2 

3 n + D n 


= F({,f) X — ^ — — cos nu (t - r ) 

n=l 


(D-87) 


The rms values of f(t) and g(t) are 
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Thus, the normalized cross-covariance becomes 
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The peak value occurs for each harmonic when r = r n . If all harmonics travel at the same 
rate, r n = constant, the maximum for the cross-covariance occurs at r n and its value is 1 . 

Thus for convective velocities that are independent of frequency, the value of the 
normalized cross-correlation is unchanged by space and frequency weighting functions. 

The point to be made is that, assuming decay in time is negligible, weighting in space 
or frequency does not change the value of the normalized covariance. Thus, a decrease in 
the normalized covariance must come as a result of. a nonlinear phase shift as a function of 
frequency for each harmonic. This may be expressed as a shifting transit time for individual 
harmonics to travel between two sensors. In the case of a nonlinear phase shift as a function 
of frequency, the peak value would be less than one, and the curve would have a greater 
width. 

Since in many constant fluid flows, the turbulent intensity does not significantly 
decrease between closely spaced sensors (i.e., the power spectrum is not a function of the 
spacing along the stream tube. 


x(t)2 = yfc E A(aj)dw >• 


the decrease and spreading of the correlation curve as a function of the separation distance 
will be a result of varying phase shifts of the harmonics and thereby varying convective 
velocities as a function of frequency. 


F. A Physical Look at Space-Time Covariance Computations 

The space-time covariance has been defined by equation (D-20) as 


c f g (r) 


= — f 

T, J 
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T, 


[f(t) - [g(t + t) - M g( - t+T )] dt 


(D-91) 


where the points of origin of the time functions f(t) and g(t + r) are spatially separated by 
some finite length % , being the time mean of the function i . 

Generally speaking, only the fluctuating component, ac, is used to preserve 
information on the structure of turbulence for correlation processes. The dc components 
only raise or lower the value of the correlation by a uniform amount. 
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A typical example is a hot wire anemometer signal. The emf representing the mean 
motion of the fluid particles at the location of the wire is the mean value of the signal. 
Equation (D-91) states that the mean value is subtracted from each representative signal. 
Thus, the space-time covariance has rejected the basic method of measuring mean flow 
velocities. 

The remaining portion of the signal, which is used for analysis, is the fluctuations 
about the mean. Probability statistics tell something about the nature of what is being 
measured with the space-time correlation. A perturbation from the mean of ±Af(t) , that 
has traversed the separation of the two probes in time, r m , with some attenuation or 

modification. A, now becomes g(t + r) = ±AAf(t) . The product of f(t) and g(t + r) 
becomes equal to AAf(t), 2 . If the perturbation is over a period of time, At, , the 
contribution to the covariance will be 


At,f(t), g(t + r), = AAf(t), 2 At, 


(D-92) 


Since the average value of f(t) is zero (the means have been subtracted out) and 
since the perturbation, Af(t), , acted over a time, At, , a cancelling oscillation 
of Af(t) 2 acting over a time At 2 must exist such that Af(t),At, + Af(t) 2 At 2 = 0 . The 
contribution to the covariance by Af(t) 2 would be AAf(t) 2 2 At 2 . Thus, AAf(t), 2 At, = 


AAf (t) 2 2 


At 2 2 
At, 


, which is the contribution to the correlation by the perturbation Af(t), , 


the contribution by Af(t) 2 being AAf(t 2 ) 2 At 2 . If Af(t), > Af(t) 2 , then At, < At 2 . 
The ratio of the contributions becomes 


AAf(t)i 2 At t _ A t 2 
AAf(t) 2 2 At 2 At, 


(D-93) 


Therefore, the larger perturbations are more heavily weighted in the covariance computation 
even though they exist for a shorter duration. 

Equation (D-92) states that the space-time covariance of temporal signals is related 
to the probability density distribution through the square of the signal from the mean. 




(a) SYMMETRICAL PROBABI LITY (b) SKEWED PROBABI LITY 
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Therefore, there is a weighting of the wings of the probability distribution. If the 
probability of. the velocity distribution is symmetrical, speeds moving faster than the mean 
velocity are weighted equally to those of a lesser velocity. The correlation would then build 
equally at the corresponding time delays, the average value being the mean velocity 
indicated by the mean value of the probability density distribution. If, however, the velocity 
distribution is skewed (Fig. (b) above), the longer wing is weighted much more heavily. The 
covariance would indicated a speed commensurate with those associated with the long wing. 

Therefore, either natural skewness of the distribution or the effects of two separate 
flows would give an indication of a velocity other than the mean. Intermittency of the 
boundary layer is just such a situation. In the outer poriton of the boundary layer, the 
model presented indicates that wave packets propagate to the outer region and are travelling 
well below the mean velocity. These wave packets are then the tails of the probability curve. 
The velocity measured from the space-time covariance calculation would indicate the 
packet’s speed rather than the mean velocity, as indicated above. In fact, there appears to be 
no reason at this time to expect the measured convective velocity to be near the mean 
velocity except at the point of production where the probability distribution will, be 
symmetrical and at the location where the turbulent packets are in equilibrium with the 
surrounding flow. 


G. Conclusions 

The values which determine the mean statistical averages of a signal are eliminated 
from covariance computations, and thus the ability to measure mean values is coincident, 
rising from the shape of the probability distribution. Eddy motion and intermittency both 
affect the shape of the distribution and thus affect the measured value of convective 
velocity. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, March 30, 1973 
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